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1. Introduction 


One of the most potential applications of solar energy is the 
supply of hot air for the drying of agricultural, textile and marine 
products, and heating of buildings to maintain a comfortable 
environment especially in the winter season. Several designs for 


Nomenclature 


area element 
height of solar air collector 
specific heat 
depth of air channel 
h equivalent diameter of air channel 
heat transfer coefficient 
k thermal conductivity 
L Length of air heater 
M mass per unit collector area 
m mass flow rate 
Nu Nusselt number 
Pr Prandtl number 
R recycle ratio 
r fraction of the mass flow rate 
Ra Rayleigh number 
Re Reynolds number 
S incident solar radiation 
T 
t 
U 
V 
W 


TUUANA ws 


temperature 

time 

loss coefficient 

wind speed 

width of absorber (collector width) 
x,y space coordinate 


Greek letters 

a solar absorbance 

ô thickness 

E emissivity 

o Stefan-Boltzmann constant 
T Transmitivity 

p reflectance 

PF fluid density 


a ambient 

b bottom plate 
c convection 

f fluid (air) 

g cover 

i inlet 

o outlet 

p absorber plate 
r radiatif 

T,t top 


solar air heaters have been proposed and discussed in literature. 
The designer and potential user of these systems must consider a 
number of factors when comparing their merits. These can mainly 
be categorized as: (i) thermal performance, (ii) cost and (iii) 
lifetime, durability, maintenance and ease of installation. Thermal 
performance of collectors is compared by using the concept of 
thermal efficiency. It is generally believed that the thermal 
efficiency of a solar collector is the major requirement for the 
prediction of thermal performance of the complete solar system of 
which the solar air collector is a part [1]. It is, therefore, important 
that this thermal efficiency information be in a form which is 
useful for existing and projected future solar energy system design 
methods. To save energy, the method of exergy analysis is useful in 
solar systems, because it appears to be an essential tool for system 
design, analysis and optimization of thermal systems [2]. One 
indication is that there are about two review papers in the area of 
solar air heaters published in recent years. The first was proposed 
by Chandra and Sodha [1]. These authors have provided a 
fundamental understanding of testing procedures for solar air 
heaters. The second was proposed by Ekechukwu and Norton [3]. 
They have classified solar air heaters broadly into two types: bare 
plate and cover plate solar energy collectors. Based on this 
classification, these authors have summarized various designs of 
solar air heaters. 

With the development of computer, hardware and numerical 
methodology, advanced mathematical models are being used to 
carry out critical investigations on solar air heaters. The advantages 
of these methods are that they can produce extremely large 
volumes of results at virtually no added expense and it is very 
cheap to perform parametric studies to optimize equipment 
performance. The second reason for such work on numerical 
simulation is that some parameters are difficult to test, and 
experimental study is expensive as well as time consuming. 

The purpose of this work is to review the present state of 
mathematical modeling of solar air heaters. The validation is an 
important step in mathematical modeling development, and 
therefore comparisons with actual experimental values or theore- 
tical results have been included where possible. Since most of the 
models have to be solved numerically, the numerical techniques 
employed in the solution of different models are presented. 


2. Mathematical models 
2.1. A bare plate solar air heater 


This consists of a single channel design with single air flow 
between absorber and bottom plates with insulation (see Fig. 1). 
Choudhury et al. [4] studied the thermal performance of this 
model. The energy balance equations for the absorber plate and the 
rear plate have been written as 


Bare plate: &œpS= hpa(Tp— Ta) + hp(Tp-— Th) + Nps (Tp — Tf) 


(1) 


Rear plate: = Mp,(Tp — Ty) = My p (Th — Tf) + Up (Tp — Ta) (2) 
With a design selection methodology whereby, for combined fixed 
values of pressure drop, air-flow rate and duct length, duct-depth 


values and the resultant values of forced convective heat-transfer 
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Fig. 1. A schematic view of a bare plate solar air heater [4,7-9]. 


coefficients, air-temperature increment per unit incident flux and 
system efficiency have been computed. From the discussion, it was 
shown that the computation methodology and the design curves 
obtained can be successfully used to construct solar air heaters 
with pre-determined acceptable values of pressure drop with the 
required rate of air flow, minimum material cost and highest 
possible efficiency. 

Thermal performance of this design has also been predicted by 
Ong [7]. The mathematical model proposed by this author is 
similar to the one of Choudhury et al. [4] with a small change in 
Eq. (1) where, hpa is replaced by Up = hw + h,gs. In this mathematical 
model, the solar collector was assumed sufficiently short for which 
the assumptions were valid. The mathematical solution procedure 
involved a matrix inversion of the mean temperature vector 
derived from energy equations. Predicted temperatures for the 
solar air heaters were presented. The effects of wind and film heat 
transfer coefficients on the prediction were discussed. 

Njomo [8] and Njomo and Daguenet [9] have investigated heat 
transfer in the above design (Fig. 1). With general assumptions in 
the heat transfer modeling, the energy balance equations of 
different components of the solar air heater have been written as: 


For selective absorber : M Cop ae = ApS — Ny po(T p — Tp) 
— hepf(Tp — T f) — (Hepa + hrpa)(Tp — Ta) (3) 


Tp | mMsCps OT py 


For air flow : bse 5Cos ae + W x 


= hepf(Tp-Tf)+hoapf(Th-Tp) (4) 
For bottom plate: MpCpp 7 
= hy pp(T p — Tp) — he p(Tp — T f) — Up (Tp — Ta) (5) 


With the mathematical model proposed, Njomo [8] predicted 
the thermal performance of such a collector. To determine a solar 
air heater thermal behavior, Njomo and Daguenet [9] used a 
sensitivity analysis which is an excellent technique to help 
designers and decision makers in anticipating and preparing for 
a solar system project. The results have demonstrated the 
existence of an important dimensionless parameter referred to 
as thermal performance factor of the collector that compares the 
useful energy that can be extracted from the heater to its overall 
heat losses. A detailed sensitivity analysis has been performed, and 
the results successfully summarized to measure the impact of 
changes in operation and meteorological parameters on the 
simulation values of the thermal efficiency and the fluid 
temperature rise between entrance and exit of the heater. 


2.2. Single one cover solar air heater 


This design consists of single channel with single air flow 
between the cover and absorber plate (see Fig. 2). Many authors have 
investigated the thermal performance of this design [7,9,13-18]. 

Garg et al. [13] developed the transient heat transfer model of a 
single one cover solar air heater. This transient model has included 


absorber plate 


Fig. 2. A schematic view of a single cover solar air heater [7,9,13,16]. 


the thermal capacitance effects and the conduction losses. The 
energy balance for each component of the collector yields the 
following: 


2 
Cover glass: gS = MgCg d + kgôg ots + hw(Tg — Ta) 
+ rgs(Tg — Ts) + Meg ¢(Tg — T f) + Mrpg(Tg — Tp) (6) 
Absorber plate: &œpTpS = M,C op kô Tp 
p : ptp PLP Be T “POP g 
+hepf(Tp -T f) + hrpg(Tp— Tg) +Ur(Tp — Ta) (7) 
. aT, GyCp Ty 
Air flow : MC-a + a z2 
= Neg p(T ¢ — Tg) + Rep f(T p — Tf) (8) 


The solar radiation data and ambient temperature have been 
represented by Fourier series as follows: 


S = 237.17 + 388.7 cos(t@, — 3.272) + 195.8 cos(2ta, 
— 0.2757) + 27.62 cos(3tm, — 3.8) + 38.44 cos(4ta, 
— 3.44) + 17.76 cos(5tw, — 4.536 (9) 


and 


Ta = 14.354 + 6.11 cos(tm, — 3.858) + 1.88 cos(2tw, 
— 0.509) + 0.225 cos(3tw; — 1.25) + 0.288 cos(4ta, 
— 2.62) + 0.146 cos(5tw, — 3.92 (10) 
with œw = (1/12) rd/h. 
The initial conditions considered by these authors were: 


T(t = 0) =T,(t = 0) =T,(t=0)=T,(t=0) and the boundary condi- 
tions were given as follow: 


OTs 7 OTs OT p -o 
OX |o OX |e ; OX |o ; 
OT, 7 
ala 0, (11) 


Eqs. (6)-(8) have been solved numerically using explicit finite 
difference analysis. The results showed the effect of various design 
parameters of the heater on its performance. 

This design of solar air heater has also been investigated by Ong 
[7]. Using the thermal network, the following equations have been 
written for a single one cover solar air heater [7]: 


AgS + rpg (T p — Tg) + Meg p(T ¢ — Tg) 
= (hw + Irs) (Tg — Ta) (12) 


Cover : 


The wind heat transfer hy, has been computed using (A1): 
Fluid between cover and absorber: ley f(T p — Tf) 


(T fm —T fi) 
= Neg p(T 5 WE 


Tg) + 2mcp (13a) 
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with 
(T jo +T 
T fin = tae (13b) 
Absorber plate: = Tg@pS = Mcp ¢(Tp — Tf) + Mrpg(T p — Tg) 
+ U,(Tp — Ta) (14) 


The methodology used for this design of solar air heater is 
similar to the one used by the authors to study a bare plate solar 
heater. 

The mathematical model proposed by Naphon and Kongtragool 
[17] is similar to the one of Ong [7] instead of Eq. (14) which is 
replaced by 


dT 
MCp— l= hye (Tp Tg) + hyp(Tp —T yf) (15) 


The above equations proposed by Naphon and Kongtragool [17] 
have been solved by the explicit method of the finite difference 
scheme. The results obtained from the model have been compared 
with the experimental data of Yeh et al. [25]. It appeared that 
numerical results and experimental data are in good agreement. 

Mohamad [15] investigated also this design of solar air heater. 
The equations proposed by this author are similar to those of 
Naphon and Kongtragool [17], instead of equation for the cover 
which did not contain the term of radiative exchange between the 
glass cover and an ambient. The main objective of this author was 
to minimize heat losses from the front cover of the collector and to 
maximize heat extraction from the absorber. 

With the aim to determine the flow channel optimum geometry 
of this single cover solar air heater (Fig. 2), Hegazy [16] proposed a 
new mathematical model. This model was based on the following 
assumptions: (i) the flow channel was hydraulically smooth; (ii) 
there was no air leakage to or from the heater; (iii) heat capacity 
effects of glass cover, enclosed air, absorber and bottom plate were 
negligible; (iv) the temperatures of glass, absorber and bottom 
plates were varied only along the x-direction of the air flow; and (v) 
the solar air heater consisted of a number elements so that the 
temperatures of the element surfaces are uniform, while that of the 
air stream inside varies linearly along its small length. With these 
assumptions, the energy balance equations have been written as 


For glass cover: gS + hrpg(Tp — Tg) + Neg f(T fm — Tz) 
= (hw + hrs)(Tg — Ta) (16) 
where Tin is given by Eq. (13b). 


For air stream between cover and absorber 


T fm —T 
Reps(T p =T fin) = Nog f(T jm — Tg) 4 26, | “Td P a7 
G is a specific mass rate. 
At absorber plate: = Tg@pS = hepf(Tp-—T fm) + hrpg(Tp — Tz) 
+ Up(Tp — Ta) (18) 


An iterative procedure has been employed by the author to 
analyze the performance of this solar air heater in terms of useful 
energy, thermal efficiency and air temperature rise. 

Hegazy [26] investigated the effect of variation in the absorber 
width on both the thermal and hydraulic performances of a single 
cover solar air heater (Fig. 2). In addition to the assumption made 
by the author [16], it was assumed that heat transfer is quasi- 
steady and one dimensional, and that turbulent air flow inside 
small element temperatures varies linearly along its small length. 
Hence, the energy balance equations of this design of solar air 


collector are given by Eqs. (16) and (18) for glass cover and 
absorber plate. For the air stream, the energy equation has written 
as 


Reps (Tp — T fm)Wdx = Neg f(T fm — Tg)Wdx + MC p(T fo — T fi) (19) 


U, has been computed using the same equation as Hegazy [16], 
in which W = W(x) is the average width of an element dx at a 
distance x from the inlet and G = m/A,. This study was considered 
as a more general and fundamental case when variable width 
absorbers were utilized, instead of conventional rectangular plates, 
in collecting incident solar radiation. The results indicated the 
variable width collectors with shape parameters exhibit perfor- 
mance behavior. 

Aboul-Enein et al. [27] have analyzed this solar air heater device 
with and without thermal storage. In order to write the energy 
balance equations, these authors have made the following 
assumptions: (i) the heat capacities of the glass cover, absorbing 
plate and insulation were negligible; (ii) there was no temperature 
gradient across the thickness of the glass cover and the storage 
material has an average temperature at time t. (iii) the system was 
perfectly insulated and there was no air leakage from the collector; 
(iv) no stratification existed perpendicular to the air flow. Without 
the storage material, the energy balance equations of the glass cover 
and for the plate are similar to those proposed in the mathematical 
model of Hegazy [16]. For the air flow, the model was written as 


aT aT 
ep s(Tp — T p)W = mC pF + djp WCs E + Inee(T y — Te) W 
+U(T p= Ta)W (20) 


with the storage material, the energy balance equation for the glass 
cover has the same form as without the storage material. However, 
the energy balance equations for the absorber plate and storage 
material were written as follow: 


TgQ@ pS = hrpg(T ps — Tgs) + hepf(Tps — T ps) 
k 

+ —? (T ps — Tst) (21) 
Zp 


For the absorber plate : 


: k 
For the storage material : ze (Tps — Tst)Ap 
p 


dT. 
= Mst Cst a t Up (Ab t As)(Tst Ta) (22) 


The forced convective heat hega is given in Appendix A. 

A computer program was prepared for the solution of this 
mathematical model. An optimization process for a flat plate solar 
air heater with and without thermal storage was carried out. 
Comparisons between experimental and theoretical results 
indicated that the proposed mathematical model could be used 
for estimating the thermal performance of flat plate solar air 
heaters with reasonable accuracy. 

Zhai et al. [18] have recently investigated a single cover solar air 
heater integrated with the building. Based on some assumptions, 
the equations governing the performance of the solar air heater 
were given by equations similar to (12)-(14) without the term of 
the radiative exchange between cover and the sky, and with 
Qi = mC,(Tyo — T,) and Up = [1/Npattic + dplkp]'. 

Njomo [14] and Njomo and Daguenet [9] analyzed the heat 
transfer in the single cover solar air heater. The following energy 
balance equations have been written for the components of this 
design: 

For combined plastic—glass transparent cover: 


dT, 
MgC gap = gS + Mrpg(T p — Tg) + Neg f(T f — Tg) 


~ (hega + hrca) (Tg ~ Ta) (23) 


1738 R. Tchinda/ Renewable and Sustainable Energy Reviews 13 (2009) 1734-1759 


oT, mC,f OT 
PreiCor ae + WOE 


=Neps(Tp —T yf) — Neg s(T ¢ — Tg) (24) 


For air flow : 


For absorber plate : Mice? = (@T)S — Nep f(T p -T f) 
hy pg(T p — Tg) — (Ub + Mr pa)(T p — Ta) (25) 


The optical efficiency of the collector is given by [29]: 


To 
(ta) = BrP 


BEETA ee) 


Tsky was given by the correlation proposed by Clark [12]. This 
mathematical model has been used by Njomo [30] to study the cost 
of the energy extracted from these collectors as a function of their 
lifetime. Numerical calculations have been performed for glass 
cover and plastic cover air collectors of the same dimensions under 
identical conditions of insolation and mass flow rate. The results 
showed that a significant gap exists between the thermal 
performances of the glass cover and plastic cover air collectors 
and the useful energy delivered by the plastic cover air collector is 
cheaper than the energy extracted from the glass cover air 
collector. 


2.3. Single two covers solar air heater 


This collector is made of two covers and a single channel air 
flow between the second cover and absorber plate, and insulation 
(Fig. 3). This model has been studied by Njomo and Daguenet [9], 
Mohamad [15], Naphon and Kongtragool [17] and Njomo [31]. 

Mohamad [15] investigated heat transfer in this design. Under 
steady state operating conditions, the energy balances equations 
for the mentioned collector written by this author are 


For the top cover: @g1S = Mw(Tgi — Ta) + hcgig2(Tg1 — Tg2) 
Tr h;g1g2(Tg1 = To2) (27) 


For second glass cover: @g2TgiS = heg2 f(Tg2 — T f) 
+ hegig2(Tg2 — Tg1) + hrgig2(Tg2 — Tg1) + hrg2p(Tg2 — Tp) (28) 


. dT 
For air stream : mCp -74 = Nego f(Tg2 -T f) + Mep (Tp — Tf) 
(29) 
For absorber plate: QpTgiTg2S 
= Neps (Tp = Tf) F h;pg2(Tp E Tg2) ar U,(T p = Ta) (30) 


Naphon and Kongtragool [17] have studied this design. The 
mathematical model obtained by these authors is similar to the 
model of Mohamad [15] with change in the equation describing 
the heat transfer in the upper cover in which the radiative 
exchange between the upper cover and ambient is taken into 
account. 


Fig. 3. A schematic view of a single cover solar air heater [7,9,16,31]. 


Njomo and Daguenet [9] and Njomo [31] have also investigated 
heat transfer in this design (with plexigas—glass cover, glass cover- 
glass cover and plastique-vitre). These authors in their investiga- 
tion have considered the two covers as a combined system which is 
characterized by one heat transfer equation. In this case the energy 
balance equations of the different component of the solar air heater 
design are given by Eqs. (23)-(25) without the radiative exchange 
between absorber plate and sky. (ta) is given by Eq. (26) with: 


Tg1Tg2 

t= 31a 

s 1- Pg1Pg2 ( 

Pz=P plati (31b) 
= s 1 — Pgi Pg2 

Og = 1 — Pg — Tg (31c) 


In the mathematical model studied by Njomo and Daguenet [9] 
and Njomo [31], the heat exchange between the two covers was 
neglected. 


2.4. Back-pass solar air heater 


In this design, the absorber plate is placed directly behind the 
transparent cover with a layer of static air separating it from the 
cover (Fig. 4). The air to be heated flows between the inner surface 
of the absorber and the layer of insulation [3]. Heat transfer in this 
design has been investigated by Choudhury and Garg [4], Ong [7], 
Hegazy [26], Garg et al. [32], Choudhury et al. [33], Al-Kamil and 
Al-Ghareeb [34], Jannot and Coulibaly [35], Hegazy [36]. 

Garg et al. [32] developed a simple mathematical model to 
investigate the performance of the system under different design 
and flow conditions. The energy balance equations for different 
components of the air heater were written as follows: 


Tg) + hrpg(Tp 


For absorber plate: (&œT)S = hepf(Tp — T fm) + hrps(Tp — Ta) 
+ h;pb(Tp — To) + hrpg(Tp — Te) (33) 


For cover:  hgs(Ts Tg) = Ut(Tg — Ta) (32) 


For the flow air: hepf(Tp—T fm) 
= hen f(T fm — Tp) + mCa(T fo — Ta) (34) 


Bottom plate: h,pp(Tp — Tp) = he f(Tp — T fm) + Ub(Tp — Ta) 


(35) 


With a deal to compare theoretical parametric analysis of solar 
air heating collector with and without packing in flow passage 
above the rear plate, Choudhury et al. [4] investigated this design 
of solar air heater (Fig. 4). The steady state energy balance 
equations for different components are presented below: 


For cover: gS + hrpg(T p — Tg) + hcpg(Tp — Tg) 
= (hw + hrs) (Tg — Ta) (36) 
I S 

h k Cover 
—_ NS EEEE S, p: 

ting K ier 
rpg E abdsor' er plate 

Flow = hepi Arep 


— o 


Fig. 4. A schematic view of a single cover solar air heater [7,16,26,31,32,36]. 
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For absorber plate: &pTgS = hep ¢(T p — Tf) + hcpg(Tp — Tg) 


+ h;pb(Tp — Ty) + Arpg(T p — Tg) (37) 

P mC p dT; 

For the flow air: hep s(Tp —T ¢) = Aen p(T ¢ — Ty) 4 W d 
(38) 

Bottom plate $ heap(Tp Tp) t hy pp (T p Tp) = Up(Tp Ta) 
(39) 


The boundary conditions for the above sets of equations were 
T(x = 0) = Tf, T(x = L) = Tyo. 

In the mathematical model written by Ong [7], there is no heat 
extracted from the system between the top cover and the stagnant 
air cavity. The air has been assumed to recirculate within the space 
under free convection in an inclined wall cavity. Hence the energy 
balance equations obtained for this design are the same as 
predicted by Choudhury and Garg [4] instead of Eq. (38) which has 
been replaced by 


2mC 
Tp) = hen p(T — To) +g TTA (40) 


Hep f(T p 


In the investigation of Choudhury et al. [4] Eqs. (36), (37) and 
(39) should be considered without the terms of natural convection 
between cover-absorber plate and radiative exchange between 
cover-ambient. Eq. (37) is replaced by 


mC f dT; 
W ar "sl. 


Tf) j hp (Tp Tf) (41) 

This solar air heater model has also been investigated 
experimentally and theoretically by Al-Kamil and Al-Ghareeb 
[34]. They studied experimentally and theoretically the effect of 
solar radiation on the solar collector. The balance energy equations 
based on the subdivision concept were developed and for all the 
components of the solar air heater, were written as 


Flowing in the channel: fey ¢(T pi — T Ai 


= hof (T fi — TriJAi + MC f(T fa-1) — T fi) (42) 

Absorber: (a@T)S+ hep ¢(T pi — Tf) + Anc(T pi — Tgi) 
Ts Ts; 

+ My pp (T pi — Tpi) + hrpg(T pi — Tgi) = PV eC oa Ar 43) 
Glass cover:  OgSFs,FgAj + hrpg(T pi — Tgi)Ai + hnc(T pi — Tgi)Ai 

= (hw + Nys)(T gi — Ta)Ai 44) 
For bottom plate: hepp(T fi — Tri)Ai + Aba(Ta — Tpi)Ai 

+ Ay pp (T pi — Tyi)Ai = 0 45) 


where F,;, and Fg are shading and dust factors. The balance energy 
equations were solved by Gauss-Seidel iterations using a 
computer program written in BASIC. The heat transfer coefficients 
were determined by familiar methods available in the literature. 
The results have shown a closing agreement between predicted 
and measured data warrant the use of the subdivision method to 
predict some useful information about the radiation effect. 
Another mathematical model for this design of solar air heater 
with a plastic film has been proposed by Jannot and Coulibaly [35]. 
They studied theoretically the evolution of the radiative fluxes 
inside the collector by a ray tracing method applied to each flux 
considering the different and successive reflections and absorp- 
tions occurring on absorber plate and on the transparent cover. The 
thermal analysis will result in the radiative balances of these two 


components. Assuming that the collector area is equal to 1 m°, the 
net radiative flux received by these two components are stated as 
For absorber: 


b = Q ps T gs i Xgi® pi 4 
P 1- pg — Ops)”  1-pgl-ap) 5 
Tgi® pi Pgi®ci 4 
H H oT 
T= Pal Ap) O Tpl Ap) P 
T= T 
o p b oa T+ 46 
(Tommy — jay 1 PT (46) 
For transparent cover: 
OgsTes(1 — Ops) A giMgi(1 — Q pi) 
Do = Stix + SBS ps) çq Osis pi) T4 
2 z 1- Pes(1 a Ops) 1- Poi (1 = Q pi) = 
Tgigi(1 — Œ pi) O pi gi 4 
H + Ogi | Psky + 
T= pg api) 8 | O T- palp) P 
— 2odlgTe (47) 


in the case of a transparent cover opaque to IR, t;=0 and 
Pgi = 1 — Og;. Eqs. (46) and (47) become: 


P, = Q psTgs Sig T- Tg 
1 — Pgs(1 — & ps) 1/Q pj; — 1/ag — 1 
Te =e 
p— tb 
"Tita l/a 1 (49) 
and 
QgsTgs(1 — Q ps) T$ -T4 
D, = Stig, + LETS B) gag p ‘sg 
7 S T= pell- Ops) 1/a pi — 1/ag — 1 
+ OQgs(EaTa — Te) (49) 


The results showed that the modifications induced by taking 
into account the partial transparency of the cover towards IR 
radiation may be important. 

Another investigation has been published by Hegazy [36]. In 
this model, Eq. (16) is reconsidered for the cover balance equation, 
with Tjm(=(Tp + Tjo)/2)) replaced by T, and the hg becomes hepg. 
Energy balance on the bottom plate was written as Eq. (35). Up is 
recalled in Table 1. For the others components of the design, energy 
balance equations have been written as 


For absorber plate: a@ptgS = Nep f(T p — T fm) 
+ (Repg + hrpg)(Tp — Tg) + hrpb(Tp — To) (50) 


For the air flow: hepf(Tp—T fm) 


Tm-Tf 
Tp) 4 26C, | ae 


= Ney f(T fin (51) 

In order to demonstrate the engineering accuracy of the 
proposed criterion, the author has determined the absorber mean 
temperature at the exit. To accommodate this necessity, energy 
balances (35), (50) and (51) have been applied progressively to 
each element and solved iteratively. It was concluded that the 
channel ratio (D/L) plays a decisive role in determining the rate of 
useful heat gain. 

To determine the flow channel optimum geometry and study 
the effect of the variation in the absorber width on both the 
thermal and hydraulic performance of this design, two new 
mathematical models have been published by Hegazy [17,26]. In 
these new models, all energy balance equations proposed by this 
author in 1999 are reconsidered instead of the equation of the air 


Table 1 
A review of heat transfer coefficient used in different models. 


References Solar air heaters systems 


Heat transfer coefficients 


Choudhury et al. [4] A bare plate solar air heater 


Njomo [8] & Njomo and Daguenet [9] 


Garg et al. [13] 
Ong [7] 


Single one cover solar air heater 


Naphon and Kongtragool [17] 


Hegazy [16] 


Aboul-Enein et al. [27] 


Njomo [14] and Njomo and Daguenet [9] 


Naphon and Kongtragool [17] 
Garg et al. [32] 


Single two covers solar air heater 
Back-pass solar air heater 


Choudhury and Garg [87] 
Hegazy [36] 
Jha et al. [43] Parallel-pass solar air heater: Type 1 


Ong [7] 


Hegazy [16,26] 


Forson et al. [3] 


Yel et al. [45] Parallel-pass solar air heater: Type 2 


The radiative and natural convective coefficients were summarized in Duffie and Beckmann [5] 

Apr and hp [6]: h= Nu..(1 + MD/L)k/D with Nux = 0.0182Re°*Pr°* and M = 14.3l0g10(N) — 7.9. N = L/D if O < L/D < 60 
and N=60 if L/D > 60. 

hcpa = 0.86 Re!/?Pr'/3 (kp /Lp), with Lp = 4WL/2W + 2L for rectangular plates and with 2 x 104 < Re <6 x 104 [10]. 
Duffie and Beckmann [5] suggested that the above equation is valid for Re > 10°. 

Rep and herp (See [11]): for the laminar flow (Re < 2100): Nu = h-Dh/K p = 4.9+ 

(0.0606(RePrDh/L)'* /1 + 0.0909(RePrDh/L)°’Pr°'”) for turbulent flow (Re > 2100) by [5]: Nu = h-Dh/K p = 0.0158Re°® 
with Dh = 4WD/2W + 2D 

hrpa = O€p(T} — T4y)/(Tp — Ta), where Tsy = Tal 0.787 + 0.0028(Ta — 273)]°7° [12] and Ta is the dew point temperature 
of the ambient air; it is related to its water vapor pressure emp (in millibars) and has been calculated from: 

emp = 1013.25/760 exp(20.519 — 5179.25/Ta), hop = O(T p + Tp) (T + Th)/(2/ep) — 1 

The important dimensions and the heat transfer parameters used were given in McAdams [19], Kreith and Kreider [20] 
The radiation coefficient between parallel plates set i-j is given as: hj = o(T; + Tj)(T? 4 T?)/(1/éi) + (1/8) —1 
The forced convection heat transfer coefficient between parallel plates is given: for laminar flow (Re < 2300) 

[21] hDh/k = 5.4 + 0.00190(RePrDh/L)!:7!/1 + 0.00563(RePrDh/L)'7!, for transition flow 

region (2300 < Re < 6000) [22]: (hd, /k) = 0.116(Re?/? — 125)Pr!/7[1 + (Dp /L)?/?](4/4y)°* and for turbulent flow 
(Re > 6000) [23] (hDh/k) = 0.018 Re?’ Pr’4, the heat transfer coefficient at the entrance region is computed with 
the same relation as hp, (see [4]), hrs = o& (Te — T$)/(Tg — Ta) 

The convective heat transfer coefficient is computed using [24]: Nu = (hDh/k) = 0.333 Re°® Pr!? with 

Dh = 4WB/2W + 2B 

hops and hegy are computed with [6]: h = k/Dh{0.0158Re98 + (0.00181Rep, + 2.92)exp(—0.0379x/Dh) } 

where Rep, = 2M/(1 + D/W) Dh is hydraulic diameter and M is a mass flow parameter. Up = [1 + 2(di, + D 

+ 5¢/W)][4in/ Kin + 1/hy]71, bin and ôg are the thickness of insulation and glass cover. 

The radiative heat transfer coefficient (see Ong [7]) 

hops and hen Were computed using [28]: Nu = hepgDh/k = 1 + 1.44[1 — 1708/R, cos( 8)]* + [1 — 1708(Sin(1.88))"°/Ra 
Cos(f)] + [(Ra cos(B)/5830)'/? — 1]* where [ ]* means that only positive values are taken into account. £ is 

a collector tilt angle. 

hpg (see Ong [7]) and hrgs = og (T? + Ty,)(Tz — Tsy) 

Negp, hrpa» Mepp (See Choudhury et al. [4]) 

hij = OE p((1 — Pg)TS — BeT4)/(1 — PpPg)(Tp — Tg) and hrpa = 8 pTir(T$ — T4y)/(Tp — Ta) 

hrgia = O(Tg1 + Ta)(T2, + T2)/((1/eg1) — 1), hrg (See [7]) with i= p and j=g 

The radiative heat transfer coefficients (see Ong [7]) 

The convective heat transfer coefficient from the upper plate to the flowing air and the pressure drop [37]: 

Nu = 2bh/k = Nuo + B(b/L), AP = [fo + y(b/L)(m?/p)(L/b)*], for Re < 2550 (laminar flow) Nuo = 5.385; 

B=0.0148Re; fo = 24/Re; y= 0.9; for 2550 < Re < 10 (transitional flow) Nu, = 4.4 x 10-*Re!?; 

B = 9.37 Re°*”!; fo = 0.0094, y = 2.92 Re~°"®; for 104 < Re < 10° (early turbulent flow) Nuo = 0.03Re°”4; 

B = 0.788 Re°*”'; fo = 0.059 Re °?; y=0.73. 

hepr and hep (see [17]) other heat transfer coefficients (see [5]) 

Up = [1 + (2/W)(Sin + D + S — 3g) ][ Sin] kin + 1/hw]', Nepg (see [27]) 

hos and hay is correlated by (see Biehl [38]): h = k/Dh{0.015Re38 + (0.00181Repp + 2.92)exp(—0.03795x/Dh)} 
The forced convection heat transfer coefficient between air and plate or between air and cover is obtained 

by using the equation [24]: h = 0.0293 Re°®k p /D. The natural convection heat transfer coefficient (see Wong [41]) 
The radiative heat transfer coefficients and the forced convection heat transfer between parallel plates for laminar 
flow, transition flow region and turbulent flow (see Ong [7]). For 9500 < Re < 22000, situations involving a 

large property variation [42] Nu = (hpn/k) = 0.027 Re®®Pr'/? (u / uy)’ 

hepfi» Negpts hcbf2, hcpf2 were calculated using the following equation: h = (k/Dpj){0.015Re?8 + (0.00181Re; 
+2.92)exp(—0.03795x/Dnj)}, Up = [1 + 2(ôin + Dı + Dz + dg/W) [Sin] Kin + 1/hy]~! W was replaced by the 

average width W(x) in the second model [26] 

Radiative heat transfer coefficients (see [7-9]), free convective heat transfer coefficient 

(see [16]), Ugj = (UA)eage/Ac = {(k/x)5j}W 

{Ugh =U;!— Ga + ira or Urm = (hw + a = (hegiga + irte 

where hegig2 = 1.25(Tfm1 — Thm2)°?°, hrpgi © 40T Fm /le7" +e — 1], hrp = 40T n2/lEp" +e! -—1] 

hrgig2» hw (see Ong [7] and (Al D); h;g2a O&g(Tg2m t Ta) (Tem t i) 

Ug! = Us! + hp + lskg! and Up ~ kslz' 


OVLL 
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Matrawy [49] 
Ammari [51] 


Wijeysundera et al. [57] 
Verma et al. [59] 


Naphon and Kongtragool [17] 
Zhai et al. [18] 
Yamali and Solmus [63] 


Ho et al. [66] 


Naphon [71] 
Chandra et al. [72] 
Ho et al. [76] 
Tchinda [78] 


Karim and Hawlader [83] 


Bashria et al. [84-86] 


Choudhury and Garg [88] 


Choudhury and Garg [88] 
El-Sebaii et al. [99] 
Ramadan et al. [96] 


Solar air heater with slats 


Two-pass solar air heaters 


Double-pass flat-plate solar 
air heater with recycle 


Double-pass flat plate solar air 
heater with longitudinal fins 
Triple pass solar air heaters 


Multi-pass flat-plate solar air 
heaters with external recycle 
Solar air heater with a 

compound parabolic concentrator 
v-Groove solar air collector 


Single pass solar air heater 
with packed, type 1 

Single pass solar air heater 
with packed, type 2 
Double-pass solar air heater 
with packed 


The forced convection heat-transfer coefficient between two flats plates: for laminar flow in a 
short conduct [21]: Nu, = D;,hj/k = 4.4 + 0.00398(0.7Re,PrDh/L)'°°/1 + 0.0114(Re;PrDh/L)'"' (j = 1, 2), 
while for turbulent flow [19,46] data:Nuj = (Djhj/k) = 0.018Re?8[1 + (De/L)*”] 

U= [1/(hne + Aric) + 1/(hw H Prad 

hrca = O&c(Te + Ts)(T2 + T2)(Tc — Ts /Te — Ta), hypg and hygq (see [7]), Natural convection (see [27]), inner 
surface convective coefficients (see [7-9]) 

Uni = [(1/Regi f) + (1/Megar)] |, Aya ohi; (T? H 17) (Ti +T;), haj = Cy (Ti + T;). 

The convective heat transfer coefficient between enclosed spaced (see [27]) 

The forced convective heat transfer coefficient (see [8,9]) 

The heat transfer coefficient (Niles et al. [24]): Nu = hD/k = 0.0333 Re®8Pr!/3. 

hgia = convection heat transfer from glazing due to wind (see Eq. (A1)) + radiation heat transfer 
coefficient from glass cover to sky referred to the ambient air temperature(see [8,9]). 

The global heat transfer coefficient from insulation plate to the attic: Up = [1/hebattic + Ab/kp] 
hrgasky (see [8,9]), Mreig2 (see [7]), hegig2 (see [28]), hegıfı (see [8,9]) for laminar flow and (see [64]): 
Nugı fı = (f 1 /8)(Re fı — 1000)Pr/1 + 12.7( f ¢,/8)°°(Pr°®” — 1) with f jı = ,/0.79 In Re fı — 1.64. 
hı = h}, hz = h}. The forced convection heat transfer between two plates (see Kays [46]), for a 
short conduit hy and hp (see [19]):Nu, p = hjDe/k = 0.0158 Rea + (De De (j = 1,2), for the 
laminar flow (see Heaton et al. [21]), hr (see [67,68] and Eq. (A4)). 

1/hg = 1/hgis + 1/hy +lg/kg 1 /Ugis = 1/Hygas + hw + 1/Hegig2 + hrgig2 OF 1/Ugis = 1/Ur + 1/h,pg1 + hı 
where Negice2 = 1.25(Tg1m — Tg2m)°?° (see [67]), hrpgi = 40TP,,/[(1/Ep) + (1/ép) — 1] 

hypr = 40T3,,/((1/€p) + (1/ér) — 1], those between two glass covers (see Ong [7]) and from 

cover 2 to ambient: hygas = &g(Tzy + Tey)(Tg2 + Tsy). 

The convective heat transfer coefficient between channel (see [7,17]). 


The radiative heat transfer coefficient between two surfaces i and j (see Ong [7]) and the 
convective heat transfer coefficient (See [73]): Nu = h,D/k = 0.036 Pr®®pr°® 
hy = h}, h2 = h}, h3 = hy, andhg = h4, . Other coefficients (see Ho et al. [66]) 


ha = 0(T2 + Tz)(Tc + Tv) /1/€y + Ay /Ac(1/éc — 1)Ac/Ap, h2 = oec (Tè + T3)(Tc + Ts)Ac/Ap, hyc = (3.25 
+0.0085T, — Tc /2D,)Ay/Ap (see [79]) with D, = 21,(e¢+ epy)/l, + ef + epv Up (see [5,11,81]) 
The convective heat transfer coefficient within the airflow duct (see [61]): Nu = Nuo + n#b/L, For laminar 


flow (Re < 2800): Nuo = 2.821, £ = 0.126 Re; For transition flow (2800 < Re < 10°): Nuo = 1.9 x 10-°Re'”9, 6 = 225; 


For turbulent flow (104 < Re < 10°): Nuo = 0.030210 °Re®’*, B = 0.242 Re°”*. the top loss coefficient U, 
(see Klein [67]) 

The radiative coefficient from absorber to the cover [52,87]: 

h, = o(T + T2)(Tp + Tg)/1/ép +1/&g —1 with &y = 2ep/1 — £p; 

The convective heat transfer (see [7]), with Dh = 2H,sin(0/2)/1 + sin(6/2) 

Hepp Men and Nepry, (ee [89]):he p f = hep p = 0.2Re0$Pr'/ °k D7, hepi f = (1/Mpry +Dp/Wk pi) |, 

with hpi = 0.255Re°$P;/*k D5! e-l, For the other heat transfer coefficients (see [5]). 

For all the heat transfer coefficients (see [88,90,91,95]) 

For the various heat transfer coefficients (see [97]) 


Negi pi NUpik fı /De = 0.2Re98Pr'/?k fı /De (see [98]), with Rep = m fDe/MAp1, De = (2/3)[eD 1 /(1 — £)], € is porosity of the 


packed bed (e < 1):€ = (Vchan — Vp1)/Vchan; Vchan and Vp; are the total volume of the channel and the volume of the 
packed bed material. Dp; = 6V’ /na (see [97]) V; = total volume of n particles selected 

randomly. Ap; = 6(1 — €)/Dp1 (see [97]). pp: = Aegis (See [43]), hepi fı = [1/h p1 f +De/D Spi] 

(see [90]) with Npı f =p pDe/k fı = (0.255/6)Re f Pr"? , Sp: is a constant and depends on the 

packed bed shape. 
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Fig. 5. (a) A schematic view of a single cover solar air heater [7,16,26,32,36]. (b) A 
schematic view of a single cover solar air heater [45]. 


flow which has been replaced respectively by Eq. (17) and: 


(T fm —T fi) 
Tp) + 2GCp Cp JAG) (52) 


with Re = 4m/P(x) and Dh = 4W(x)D/P(x). 


hp(Tp — T fin) = hep f(T fm 


2.5. Parallel-pass solar air heater 


2.5.1. Type 1 

The solar air heater configuration under consideration is 
shown in Fig. 5a. It consists of one cover, double channel 
design with double air flows between cover and absorber 
plate and between absorber and bottom plates and with 
insulation provided. Many papers have investigated this design 
[7,16,17,26,39,40,43]. 

Jha et al. [43] obtained the energy balance equations in partial 
differential forms, which govern the behavior of the system, by 
considering energy conversation at each component of the system 
separately. For the solar heater shown in Fig. 5a, the following 
equations were obtained: 


At the cover : MgCy SE = oS +t hrpg(Tp — Tg) + Ne pig(T fı — Tg) 
Hegw(Tg — Tw) — hrga(Tg — Ta) (53) 


For air flowing between cover and plate absorber : 
OT f -Cp OT py 


MaCs- wW hepfi(Tp — T f1) — hefig(T fı — Tg) 
(54) 
ar PT 
For Plate 1 : M: g = Æ pTgS — kpôp Pri hrpg(Tp — Tg) 
_ Nepsa(T p T f2) hy pp (T p Tp) hepfi (Tp — Tf) (55) 


For air flowing between absorber plate and bottom plate : 
OT f2 _ —G2Cf OT f2 


k Nep f2(T p — T f2) — ho f2(Tp — T f2) 


M f2C ¢ 


ot Ww 
(56) 
2 
For bottom plate: M,C, z kpôp ae + hy pp(T p — Tp) 
— Neb f2(Tp — T f2) — Np (Tp — Tr) (57) 


The boundary and initial conditions are given by the following 
equations: 


arp} ap} ITa 
OX |y—9 0, OX |y 0, OX |x- 9 (>8a) 


Ong [7] obtained from the thermal network the following heat 
balance equations: 


For cover: gS + hrpg(Tp — Tg) + Neg fi (T fı — Tg) 
= (hw + hrs)(Tg — Ta) (59) 


For air flow between cover and the absorber plate : 


hepfı(Tp = Tf) = Neg p(T pi 


Tp —Tj 

Tg) +2m,Cr aL i (60) 

For absorber plate: = a@ptgS = hepf2(Tp — T f2) + hrpg(Tp — Tg) 
+h,p(Tp — To) + Mep pi (Tp — Tp) (61) 


For air flow between absorber plate and the bottom plate : 


Tf- Tifa 
WL 


Nep f(T p — T f2) = ha f2(T f2 — Th) + 2M2C f (62) 
With the aim to determine the channel optimum geometry of the 
flat plate and to study the effect of the variation in the absorber 
width on both the thermal and hydraulic performances of this solar 
air heater (Fig. 5a), Hegazy [16,26] has written two mathematical 
models. In the two models, the energy balance equations on the 


glass cover, the absorber plate and the bottom plate are given as 


gS + hg(T fmt — Tg) + hrpg(Tp — Tg) = (hw + hrgs)(Tg — Ta) (63) 


O@ pT gS = Nep p(T p — T fm) + Arpg(T p — Tg) + ep fo(T p — T fm) 
+ hrpg(Tp — Tp) (64) 


hes f2(T fim2 — Th) + Mrpp(T p — Ty) = Ub(Tp — Ta) (65) 
The energy balance equation of the air flow in the upper and the 


lower channel are 
For the first model: Mp1 (Tp — T fm) 


T fm — Tifi 


= he(T fm — Ta) + 261C "Ty (66) 
T m2 — Tif2 

hpf2(Tp — T fina) = hp f2 (T fm2 — Tp) + 2G2Cp (L/dx) (67) 

For the second model: Np py(Tp — T fm) = Ag (T fm — Tg) 
T —T; 

p IGC 2 i 68 

> (A/A) FO e 

hp f(T p — T fn) = ho f2(T fa — To) + 2G2Cp A (69) 


(Ap/A(x)) 


with Re = 4m/jzP and Dh = 4W/P for model 1 [16], Re = 4m/P,(x) 
and Dp = 4W(x)/P;(x) for the second model [26]. 

To study the performance of this solar design (Fig. 5a), Naphon 
and Kongtragool [17] proposed another mathematical model. This 
model is characterized by the following equations: 


Og1S = hw(Tg1ı — Ta) + he figi (Te1 — T f1) 
+ hrgı p(Tg1 = Tp) + hragı(Tg1 — Ta) (70) 


For the cover : 


dT p 
dx 
= hg p1(Tg1 — T f1) +hepfilTp-Tf) (71) 


For the air flow between upper cover and absorber: mC, 
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For the absorber plate : œ pTg1S 
= he fip(T p — T f1) + Me gap(T p — T f2) + rei p(T p — Tg1) 


+ hy p(T p — Tp) ä 72) 
For the air flow in lower channel: mC, 
= he sop(Tp — T f2) + Repa(Tp — T f2) 73) 


For the bottom plate: hef2b(Tp — T f2) + Arpp(Tp — Tp) 
+ Up(Tp — Ta) = 0 74) 


Forson et al. [43] have investigated theoretically and experi- 
mentally heat transfer in the design shown in Fig. 5a. Making a 
general assumption, Forson et al. [43] formulated the equations 
governing the performance of the system by coupling the energy 
equations of the components of the solar air heater with those for 
the useful heat extracted in the two channels. The energy balance 
equations for various components of the system are given below: 


For the cover: Sg = hw(Tg1 — Ta) + res(Tg — Ts) 
+ Mrpg(Tg — Tp) + Neg 1 (Tg — T f1) (75) 


For the air flow between upper cover and absorber: 


dT 
mC 7 = Whey p(T p — T p1) + Wheg fi (Te — T 1) 


+ Uegsı (Ta — T f1) (76) 


For the absorber plate: 


Sp = hcpf2(Tp i T f2) +F Rep fi (Tp = Tf) F h;pb(Tp = Tp) 
+ hrpg(Tp — Tg) (77) 


For the air flow in lower channel: 


dT f2 
dx 


m2C f = Whe f2b(Tp — T f2) + Whep f2 (Tp — T f2) 
+ Ues2(Ta = T f2) (78) 
For the bottom plate: 


T f2) + Depo (Tp 


he f2b(Tp Tp) + Up(Tp — Ta) = 0 (79) 


The pressure drop is given by Achenbach and Cole [44]: 


D: = 233.9hp, (h) (80a) 
T m 
with 
5 
T fm = g [Ta — 0.65(Ta — T fo)] (80b) 
and 
h = Lsin(B) (80c) 


The mass flow rate of air in each of the two channels is given by 


2Dhp; TT. 
m; = Wo SjCyi,| ot ¢ 233.9hp, | —@ —* 81 
E na e(Te) i 


Dh; is the hydraulic mean depth of the duct and CG, the 
coefficient of the velocity. The edge heat losses are estimated by 
assuming one dimensional sideways heat flow around the 
perimeter of the collector system. In Eqs. (75) and (77) S and 


Sp are given by the following: 


Š, = aigRKy cone (82) 
Sp = 08 (+ e D) aaah (83) 
with 
Ho = A a Gsc (: + 0.033Cos e3) cos( in @)sinos 

+ aa 4 sinosn)| (84) 


A FORTRAN computer program was written and used to solve 
the above set of equations. In the solution procedure it was 
assumed that the air heater is subject to a constant insulation and 
hence, the plate temperature is constant. It was shown that 
significant improvement in solar air heater performance can be 
obtained with an appropriate choice of the collector parameters 
and the top to bottom channel depth ratio of the two ducts. The 
mass flow rate was shown to be the dominant factor in 
determining the overall efficiency of the heater. 


2.5.2. Type 2 

The energy flow diagram of such a device is presented in Fig. 5b. 
At seen in Fig. 5b, two air streams of different flow rates but with 
total flow rate fixed, are flowing steadily and simultaneously 
through two separate channels (above and below the absorbing 
plate). With the presence of the second cover, this device is 
different to the one of Fig. 5a. Yel et al. [45] investigated heat 
transfer in this design. In the mathematical model proposed by 
these authors, the following assumptions have been considered: 
the temperatures of absorbing plate, bottom plate and bulk fluids 
were functions of the flow direction only, and both the covers and 
fluid do not absorb radiant energy. The energy balances within the 
differential length dx, were given as 


For inner cover: hrpa(Tp — Tgi) + Negi pi (T fi — Tgi) 
= g1al(Tgı = Ta) 85) 


Absorbing plate: a@ptzS = Ur(Tp — Ta) + Ug(Tp — Ta) 


+ Acpr1(T p — T f1) + Acpr2(T p — T f2) 86) 
For the bottom plate: he»f2(T f2 — Th) + Mrpp(T p — Tp) 

= Up(Tp — Ta) 87) 
For fluid 1: rmCy a = Whep pi (T p — T fi) 

+ Wheg 1 (Tg — T p1) + Me pigi (T fi — Tgi) 88) 

dT 

For fluid 2: (1 —r)mCy a 

= Wherr2(T p — T f2) + Who fo (Tp — T f2) (89) 


2.6. Solar air heater with slats 


Yeh et al. [47] have presented a mathematical model for double 
flow solar air heaters with fins attached. Fig. 6a shows the design of 
a solar air heater with fins attached in which the absorbing plate 
divides the air conduit into two parts, channel 1 and 2. The energy 
flow diagram of such a device is also presented. As seen, two air 
flow of different mass flow rate but with total rate fixed, flow 
steadily and simultaneously through two such separated channels 
of the same size for heating. The method for theoretical prediction 
of collector efficiencies as well as the experimental procedure was 
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Absorber plate 


Insulation 


Fig. 6. (a) Schematic view studied by Yeh et al. [47]. (b) Schematic view studied by 
Ammari [51]. (c) Schematic view studied by Kabeel and Mecarik [50]. 


similar to that presented by the same authors [45], except that the 
allowance was for the fins to be attached to the upper and lower 
surfaces of the absorbing plate. It was shown that providing fins 
attached on the collector, will improve the collector efficiency and 
that the application of the concept of double flow in the design of a 
solar air heater with fins attached was technically and economic- 
ally feasible. 

Al-Nimr and Damseh [48], Matrawy [49], Kabeel and Mecarik 
[50] and Ammari [51] have investigated a solar air heater system 
that has metal slats connecting the absorber plate to the bottom 
plate as a development to enhance the thermal performance of the 
solar heater (Fig. 6b). 

With the aim to study the dynamic performance of this design 
of solar air heater, Al-Nimr and Damseh [48] developed a 
mathematical model on the base of the following assumptions: 
(i) there are two separate phases, fluid and solid, which are not in 
thermal equilibrium. This was the case for low thermal con- 
ductivity fluid and low volumetric heat transfer coefficient 
between the fluid and solid domains; (ii) the fluid was one 
dimensional plug flow; (iii) the fluid was laminar and Newtonian; 
(iv) the properties of the fluid and solid domains were constant and 
uniform; (v) thermal and mass diffusion were assumed to be 
negligible; (vi) transient response resulted from sudden change in 
the intensity of the incident solar radiation and/or inlet fluid 
temperature. This model is described by the following equations: 


Ot Tx. PsAsCsm f 


WLp pA fUu 
PsAsCsm f 


bLp pA fU1 
PsAsCsm f 


(As -0f) (90) 


80; df dlU, g 
ət OX 


_ BLU; 
Csm ¢ fs Csm ¢ 


(9s — 95) (91) 


where 
(Ts —Tx) (Ts — Tx) (T ¢(t,0) — Tx.) 

Of = = , O; = i Tx Ao = f Ta ; 
tmp 

T= 92 


with the following initial and boundary conditions: 0{0,x)= 
0(0,x) = 0 and 07,0) = O0(T). 

The authors have outlined analytical solutions for the fluid and 
solid domains. The effects of different design parameters on the 
thermal performance of the heater were analyzed and the validity 
of the theoretical model was verified experimentally. It was found 
that both theoretical and experimental results were in good 
agreement. 

Matrawy [49] proposed another mathematical model of this 
solar air heater design. The model was based on the following 
assumptions: (i) uniform temperature distribution across the 
width of absorber plate and of back plate; (ii) uniform distribution 
of the air flow along the collector channels; (iii) the physical 
properties of the flowing air was independent of temperature and 
pressure; (iv) only steady states were considered and the energy 
stored in the components was neglected; (v) all convective heat 
transfer coefficients between the collector components and 
flowing air were equal and constant through out the solution 
procedure (vi) the temperature of the flowing air was uniform over 
the collector cross-sectional area. The energy balance equations for 
the collector components were written as follows: 


For the absorber: U;(Tp — Ta) + Neps (Tp — Tf) + Arpp(T p — To) 
A (tani) =s f a 
For the plate h,p5(Tp — Tp) al Anks D 
= hea ¢(Th — T f) + Ub(Tp — Ta) (94) 


dT. dT 
For the metal slats : ( Ak ‘) ( Ak. ‘) 
5 dz z=0 5 dz z=D 
z=D 
= J haslet; -T dz (95) 
z=0 


For the airflow: hep s(Tp —T 5) + hep f(Tp — Tf) 


z=D 
1 
+ J 2hamkeTs —T,)dz=4q, (96) 
z=0 


where A< is the collector total surface area (WLc) and A is the slat 
cross-sectional area (ôLc). Nis the number of metal vanes per meter 
of width. The determination of the radiative heat transfer was 
given in Duffie and Beckmann [52]. The convective heat transfer 
coefficients were computed on the basis of empirical correlations 
developed by Heaton et al. [21] and Nusselt [53] and given by Ong 
[7]. The solution of the above equations was carried out using an 
iterative procedure. It was concluded that the air collector with a 
box-type absorber can easily be modeled mathematically through 
the development of expressions for its parameters. The maximum 
number of metal vanes and depth of the air duct in a solar air 
collector are interesting factors for collector design. 

The mathematical model proposed by Ammari [51] was similar 
to the above for the energy balance of the absorber plate, back 
(bottom) plate, slats metal (vanes) and air flow. For the cover 
Ammari [51] has proposed the follow: 


Tg) + hrpg(Tp 


The sky temperature is defined as [54]: Tsky = 0.05527}. The 
flow is assumed to be hydro dynamically developed at the collector 
inlet. 


Ane(T p Tg) = hw(Tg — Ta) + hrga(Tg — Ta) (97) 
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Kabeel and Mecarik [50] have investigated a longitudinal fins 
collector with two covers. This model is show in Fig. 6c. By 
considering the energy balance for each part of the system in 
steady state heat flow were written as follow: 

For the top glass: — a@g1S + Mygig2(Tg2 — Tg1) + Aegiat (Tai — Tg1) 

—hw(Tgi — Ta) = 0 (98) 


For the air between the two glasses:  heg2a1(Tg2 — Ta1) 
T Acgiai (Tg1 = Ta1) =0 (99) 


For the second glass cover: —Tg10tg2S + Mega f(T f — Tg2) 
T Aygo p(T p = Tg2) = Nygig2(Tg2 = Tg1) = heg2a1 (Tg2 = Ta) =0 
(100) 


For absorber plate : 
— hrpg2(T p — Tg2) 


Tg1 Tg2M pS = Zhp f(T p = Tp) 
Up(Tp— Ta) =0 


(101) 


Z is the ratio between the actual absorber area and the absorber 
area normal to the solar radiation falling upon the collector. 


2.7. Two pass solar air heaters 


The concept of a two-pass air heater was introduced by 
Satcunanathan and Deonarine [55] and later considered by Caouris 
et al. [56] for liquid based systems. These authors conducted 
experiments on the two-pass design with air flow between the two 
glass covers and under the absorber. Wijeysundera et al. [57] 
motivated by a need to consider in greater detail both analytically 
and experimentally the concept of two-pass design, developed 
heat transfer models for two-pass flow arrangements. In view of 
uncertainties which exist in the heat transfer coefficients for 
various processes, they are made the following assumptions to 
keep the analysis simpler: (i) The temperature, of the covers and 
absorber plate were uniform. (ii) For the thermally developing flow 
in a heater channel, suitably averaged heat transfer coefficients 
available in the literature were used. (iii) The bulk-mean 
temperature of the cooling fluid was changed only in the flow 
direction. (iv) The area of the collector was large compared to the 
thickness so that end losses and shading can be neglected. (v) 
Effects due to air leakage and temperature drop in the plenums 
were negligible. (vi) The absorber plate and the covers were 
diffuse-grey for long wave radiation. For incoming solar radiation 
the covers behave like specular reflectors. (vii) The collector was in 
quasi-steady state. Fig. 7 shows a schematic view of the two-pass 
flow arrangement under consideration. The air flows first between 
the covers and then through the bottom channel. There is a 
stagnant air gap between the absorber plate and the cover 2. For 
the various elements of the collector, Wijeysundera et al. [57] 
wrote the nodel energy equations as 


For cover 1: Qf =Acà1S + Achrgagi (Tg2 — Tg1) — Ac (hg1a 
+ hw)(Tgı — Ta) (102) 
h, 

ae Tatha NAL Cover? 

Diar N 

— hri heoi s 
aa T anana NEA a Cover] 
Tp hiper Nepet Absorber plate Tim 
f 
Tiog—— hipo jä 


Insulation 


Fig. 7. Schematic view studied by Wijeysundera et al. [57]. 


For cover 2: Q f2 = Acà1S — Achrgig2(Tg2 — Tg1) 


+ Ac (h; pg2 = Nepg2)(T p Tg2) 103) 


For the absorber plate: Qpf 
= Acà pS — Achy pp (T p —T;) — Ac(Mrpg2 + hcpg2)(T p —Tg2) 104) 
For the bottom plate: Qj) - =Achypp(Tp — Ty) — Acha f(T — Ta) 
105) 
For the air flow between covers: Qf 


1 yh 
= Aan |(T2— Tr) + F(T -T p) + FRE (Ta ~T)| 106) 


For the airflow between absorber plate and bottom plate : 


h 
Q fo = AcUa1 | (Ty — Ty) +F (Ta — Tp) + FS (Typ — T a) 
Neg f 
(107) 


with F' = mC/(Nega f + hegi f)WL[1 — exp(—(hegs f + hegi pf) WL/MC)). 
Q is the total heat flux and the fractions A; of the solar radiation 
absorbed in the various elements of the collector are calculated 
using the net method [58]. A similar procedure has also been 
applied by Wijeysundera et al. [57] on the two-pass solar air heater 
with two covers. Wijeysundera et al. [57] used this mathematical 
model to study the thermal performance of the two-pass solar air 
heater design. The results show that the performance of the two- 
pass design with air flow above and below the absorber plate in 
turn is better than the single pass system. The results of the 
experiments and predictions of this heat transfer model have 
shown that there is some scatter in the experimental data. There 
results were also compared to the experimental data of 
Satcunanathan and Deonarine [55]. 

A new mathematical model has been proposed by Verma et al. 
[59]. They wrote down the energy balance equations for the design 
of Fig. 8 under usual assumptions [13] to show the optimization 
procedure used. The two-pass solar air heater design studied by 
Verma et al. [59] is different from the one modeled by 
Wijeysundera et al. [57] (see Figs. 7 and 8). These analyses give 
the following equations: 


At cover glazing: agS = hrgs(Tg — Ts) + heg pi (Tg — T f1) 


+ hw(Tg — Ta) + hrgp(Tg — Tp) (108) 


At a fluid stream between cover and absorber : 


mC dT f 
=hepip(T fi — Tp) +94 ie 


Neg pi (Tg i Tf) 


(109) 


At absorber plate: &pTgS + hrgp(Tg — Tp) + hefip(T fi — Tp) 
= hepf2(Tp — T f2) + hrpo(Tp — To) (110) 


Cover glazing 


Bal I degr 


Nrpg 
—* S heptt Absorber Ea) 
h rpb’ i hepr 


Ncbi Bottom plate 


EE | mee _ __ _ 


Up 


Fig. 8. Schematic view studied by Verma et al. [59]. 
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At a fluid between absorber plate and a rear plate : 


Rep s2(Tp — T f2) = Ne pao (T f2 — To) + -7 d (111) 


At rear plate : 


Ty) + hrpb(Tp — Tp) = Ub(Tp — Ta) 


(112) 


Ne fap (T f2 


The differential equations can be solved by using the following 
boundary conditions: 
Trilol an and Ty.) =T paler (113) 

The temperature of the sky T; have been calculated by using the 
expressions as given Swinbank [54]. An iterative procedure has been 
used by Verma et al. [59] to solve the above Eqs. (109)-(113). In this 
procedure, the temperatures of the cover, plate, air, etc. were initially 
assumed. The new temperatures were obtained by using the 
solution equations and were compared. The procedure was repeated 
until the solution is converged and finally the collector outlet air 
temperature was obtained. It was shown that there exists an 
optimum mass flow rate corresponding to an optimum flow depth. 

The mathematical model proposed by Wijeysundera et al. [57] 
was used by Naphon and Kongtragool [17] to study numerically the 
effect of air mass flow rate on the heat transfer characteristic and 
performance of the single two-pass solar air heater. The explicit 
method of finite difference scheme was used to solve these models. 

In order to study the effects of various design and operational 
parameters (duct length L, depths of inlet D4, and outlet D> air flow 
channels and air mass flow rate) of the system, Choudhury et al. [33] 
modified a mathematical model proposed by Wijeysundera et al. 
[57]. The schematic view is identical to the one proposed by Verma 
et al. [59] (see Fig. 8). The steady-state energy balance equations for 
different plates and flowing air in the inlet and outlet channels of the 
air heaters, assuming the air inlet temperature to be the same as the 
ambient temperature, were written as follows [33]: 

For the two-pass, single-cover air heater: 


QgS = hy(Tg — Ta) + Neg f1 (Tg — T f1) 
+ higp(Tg — Tp) 114 
Fluid between cover and absorber plate: = mC(T f — Ta) 


= Neg pi (Tg -T f1) + Nep pi (Tp —T pr) 115 


Cover glazing : 


Absorber plate: &pTgS = Neppi(Tp — T f1) + hep fo(T p -T f2 
+ hy pp(T p — Tp) + Mr pg(T p — Tg) 116 


At a fluid between absorber plate and a rear plate : 
MC(To — T fo) = Meg f2(Tp — T f2) + heb f2(Tp — T f2) 117 


At rear plate: hy py(T p — Ty) — New fa(T p — T f2) 
= Uy (Tp — Ta) 118 


From the results discussed by these authors, it was concluded 
that the optimum values of the duct depths, which correspond to 
the minimum annual cost per unit solar energy gain, were different 
for different duct lengths and air mass flow rates. 

Chauhan et al. [60] studying the drying characteristics of 
coriander in a stationary 0.5 tone/batch capacity deep-bed dryer 
coupled to a double-pass solar single air heater and a rock bed 
storage unit, proposed a new mathematical model for this air 
heater design (Fig. 8). The energy balance equations for different 
components of the solar air heater were given by the following: 


At cover glazing : 


hw(Tg — Ta) 


dT, 
MgCg — E = agS — Neg fi (Tg — T f1) 


dt 
hrgp (Tg — Tp) (119) 


At a fluid stream between cover and absorber : 
mC p dT dT 
f Pg fl 
w dx MICI ae 


Negi (Tg — T f1) 


(120) 


—hepip(T fi —Tp) = 


At absorber plate: @ptgS+ hrgp(Tg — Tp) + Ae pi p(T pi — Tp) 


dT, 


— hepf2(Tp — T f2) — hrp(Tp — Tp) = MpCp -i (121) 
At a fluid between absorber plate and a rear plate : 
hepf2(Tp — T f2) — he f2(T f2 — Ty) = ri a + M Cp oe 
(122) 
At rear plate hef2b(T f2 — Ty) + Ay p(T p — Tp) — Up(Tp — Ta) 
= M,C, Ët (123) 


The initial and boundary conditions were: Tp(x,0) =T,(1), 
Tpo(x,0) = Ta(1), Tp(x,0) = T,(x,0) = T.(x,0) = Ta(1), Tp (0,t) = 
Tn (Lit) = To(0,t), To(Lt) = Tjo. 

Choudhury at al. [33] have investigated a double pass solar air 
heater with double cover. This model is shown in Fig. 9. It consists 
of double covers, double channel design with double air flows 
between inner cover and absorber plate and between absorber and 
bottom plates and with insulation provided. The mathematical 
model proposed by these authors is as follows: 


a= Tf, 


Cover 2: Og2S = Aw (Tg2 _ Ta) + heg2g1 (Tg2 — Tg1) 
ZA hrg2g1 (Tg2 a Tg1) 124 
Cover 1: Œg1 Tg2S = h;g2g1 (Ta — Tg2) + Aegag1 (Tg1 — To2) 


+ hogz fı (Tg1 —T f) + hr pgi (Tg —Tp) 125 


Fluid between absorber plate and a cover 2: mC(T fio — Ta) 
= Negi fi (Tg1 — T f1) +hepfi(Tp -T f1) 126 


Absorber plate: &ÆpTg2Tg1S = hrpgı (Tp — Tgi) + hepfi(Tp -T f) 
+ hepf2(Tp — T f2) + hrø(Tp — To) 127 


At a fluid between absorber plate and a rear plate : 


MC(T p29 — Tfio) = hegf2(Tp — Tf2) + Mebp2(Tp — Ty2) (128) 
At rear plate hy pp (T p20 _ Tp) = hef2(Tp = Ty2) 
= U, (Tp — Ta) (129) 


Choudhury et al. [33] has not taken account of the heat transfer 
by the radiation from the cover to the sky. In the above equations 
Tp, and Tp are given by: Tp, = (Tro + Ta)/2 and Tp = (Tpo + Tfo)/2. The 
various heat transfer coefficients used in this analysis have been 
computed by using the equations given in Duffie and Beckann [5], 
Tan and Charters [6] and Hollands et al. [61] (see Table 1). 
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Fig. 9. Schematic view of energy-flow of double-pass type solar air heater (see 
Choudhury et al. [33]). 
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Recently, to evaluate the effects of two kinds of roof solar 
collector, Zhai et al. [18] used the mathematical model developed 
by Choudhury et al. [33]. In this work, for the natural ventilation 
mode, the mean air temperatures in each channel were calculated 
according to Hirunlabh et al. [62]. 
T fm = 0.75 T29 + 0.25 T, (130) 


and the air mass flow rate per unit of the system was given by 


nape ees =T,) (131) 


Api T, (Ein + Eo F f5) 


with £in = 1.5, Eo = 1.0, f= 0.056. The results show that the double 
pass roof solar collector is more efficient. 

Recently Yamali and Solmus [63] have investigated a desalina- 
tion system which has a double-pass flat plate solar air heater with 
two covers as a component. In their work, a transient analysis of 
double-pass solar air heater has been carried out. The energy 
balance equations are follows: 


dT 

For cover 2: MgCg n = AgSWL + WL h,g1g2(Tg1 — Tg2) 

— WL hw(Tg2 = Ta) — WL Nygasky(Tg2 — T sky) + WL hegig2(Tg1 = Tg2) 

(132a) 
with 
Tsy =Ta — 6 (132b) 
dT 

For cover 1 : MCs -gt = AgTgSWL, + WL hy poi (Tp — Tg1) 

— WL legs 61 (Tg1 — T f1) — WL hegig2(Tg1 — Tg2) 

— WL Negig2(Tg1 = Tg2) (133) 

dT 

First air pass : MCE = WLhepfi(Tp -T fi) 

+ WLhegi f1 (T gi — T f1) — MC f(T fio — T pri) (134) 
Absorber plate : c ua 2SWL, — WLh 

sorber plate: Mg Sap T ete Le — cp f2(T p — T f2) 


= WLhep fi (Tp = Tf) = WLH; poi (T p = Tg1) = WLh,pb(Tp = Tp) 


(135) 
. dT f2 
Second air pass: MpCp dt WLhepf2(Tp — T f2) 
+ WLI ep f2(Tp — T f2) — MFC El p20 — T fio) (136) 


The four order Runge-Kutta method was used to solve 
numerically the energy balance equations. For this purpose, a 
computer simulation program based on mathematical model was 
developed by means of MATLAB software. Theoretical results 
illustrated that the system productivity is strongly affected by the 
solar air heater area and slightly influenced by the wind speed 
variations and bottom heat loss coefficients of the solar air heater. 

With the aim to reduce losses from the front cover of the 
collector, Mohamad [15] modified the two-pass solar air heater 
with two cover proposed by Choudhury et al. [33] and 
Wijeysundera et al. [57]. The schematic view is shown in 
Fig. 10. This model consists of double covers, double channel 
design with double air flows between upper cover and inner cover 
and between inner cover and absorber plate, the lower surface of 
the absorber is in contact of the insulation provided. Under steady 
state operating conditions, he has written the energy balances for 
this design as follows: 


QgoS = hw(Tg2 — Ta) + Mega pi (Tg1 — T f1) 
+ hrgig2(Tg2 — Tg) 


Cover 2 : 
(137) 


hrg g2 = 
e en A NPs cecsorecassvseesiacel C anlan) 
ha Rov 


Fig. 10. Schematic view studied by Mohamad [15]. 


Fluid between the covers: Meg ¢,(Tg2 — T f1) + Me pigi (T fi — Tgi) 


— W dx aa 

Cover 1: Ogi Tg2S = hrg2g1 (Ta = Tg2) + Negi fi (Tg1 = Tf) 
+ he fag1 (Tg1 = T f2) + higi p(T 1 = Tp) (139) 
Fluid between the cover 2 and the absorber plate : 
mC s dT 

he ogi (Tet — T y2) + Resap(Tp -T = GE (140) 
For the absorber plate: Q@pTg2TgiS 

= Ne pop (Tp — T f2) + hrpgı(Tp — Tg1) + Up(T p — Ta) (141) 


In order to increase heat transfer from the absorber to the air, 
Mohamad [15] suggested replacing the absorber plate by a 
porous matrix absorber. In this case, if it was assumed that the 
air and the solid matrix were in thermal equilibrium (the 
temperature of solid was equal to the temperature of the air 
locally), the energy balance in the second passage can be 
simplified as 


2 
T p2) + keff + pT et TeS 


h f2g1(Tg1 — T f2) + Ua(Ta 


mC p dT 
= fe fe 
=- W dx (142) 

The results indicate that, under normal operating conditions, 
the thermal efficiency of the suggested heater is much higher than 
the efficiency of the conventional air heaters. 


2.8. Double-pass flat-plate solar air heater with recycle 


Recently, it was reported that increasing the fluid velocity by 
using recycle in double-pass parallel-plate heat exchangers 
enhances the heat-transfer coefficient, resulting in improved 
performance [65]. The schematic view of the design of the solar 
heater is shown in Fig. 11. In addition to increasing the fluid 
velocity, the recycle operation also produces an effect by 
remixing the inlet fluid with the hot outgoing fluid. Due to 
the fact that the strength of forced heat convection by using 
recycle effects plays an important role in designing and 
operating the solar air heaters, Ho et al. [66] developed the 
mathematical formulation for a double-pass solar air heater with 
recycle. In this model, the following assumptions were made: the 
temperatures of the absorbing plate, bottom plate and bulk fluid 
are functions of the flow direction only. Both glass covers and 
fluids do not absorb radiant energy. The radiant energy absorbed 
by the outlet cover is assumed to be negligible [66]. With the 
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Fig. 11. Schematic view of energy-flow of double-pass type solar air heater with 
recycle (see Ho et al. [66]). 


above assumptions, the energy equations were written as 
follows: 


Cover 1: Mrpg2(T p — Tg1) + hy (Tp(Z) — Tg1) = Ugis(Tg1 — Ts) 
143 
Absorbing plate: &ŒpTg1ıTg2S — hr(Tp — Ts) + hg(T p — Ts) 
+ hy (Tp — Ty(x)) + ho(Tp — Ta(X)) 144 
Bottom plate: h,pr(Tp — Tr) + h, (Ta(x) — Tr) 
= Upis(Tp — Ts) 1i 
Channela: hz(Tp-— Ta(x)) — hy (Ta(x) — Tr) 
(R-+1)mCp] dT q(x) 
= | a 5 146 
Channel b: hı(Tp — Ty(x)) — hi (Ty (x) — Tg1) 
_ [—(R+ 1)mC,] dTp(x) 
| y | Hl (147) 


In this equation, the characteristic length is the equivalent 
diameter of the duct. For the open conduit, lower channel and 
upper channel, they are respectively given by: 


_ 4HW _ (4AH)W 
Deo = ZHL W’ ea = QAH) +W and 
_ (40.—A)H\w 
Dea = (200 — A)H] +W (148) 


and the air velocities of the open conduit, lower channel and upper 
channel used by Ho et al. [66] are given by 
E m - (1+R)m 


: : [m(1 + R)] 
o = QHW)p’ a= (2AH)Wp 


RC- A)H|Wp 
(149) 


and = 


A is the ration of channel thickness. 

The experiment had been carried out by these authors. With an 
iterative procedure, Ho et al. [66] investigated the effect on the 
collector efficiency with ratio as the variable parameter for this 
type of solar air heater design (see Fig. 11). The results showed that 
the theoretical prediction agree reasonably well with experimental 
data. It was also confirmed that the collector modules with 
external recycle are thermo-hydraulically better. This is economic- 
ally feasible in the design of double-pass flat solar air heaters. 


2.9. Double-pass flat plate solar air heater with porous media 


Naphon [69] studied theoretically the heat transfer character- 
istics and performance of double-pass flat solar air heater with and 
without porous media. The model was based on that of Naphon 
and Kongtragool [17] with the following assumption: (i) flow of air 
is steady, (ii) outside convective heat transfer coefficient is 
constant along the length of solar air heater, (iii) inside convective 


heat transfer coefficient is constant along the length of solar air 
heater, (iv) thermal conductivity of porous media is constant along 
the length of solar air heater. 

The mathematical model proposed in this study for the design of 
the double-pass single solar air collector with porous media is based 
on the one proposed by Naphon and Kongtragool [17] and 
Wijeysundera et al. [57]. The equations of the top glass cover, the 
first-pass air, the absorber plate and the bottom plate are similar to 
given by Eqs. (108), (110) and (112). With the presence of the porous 
media in the second-pass air, heat transfer Eq. (111) becomes: 


dx 


aT 
ra + he f2p(T p — T f2) + Ne fan(To — T f2) 


mC p =Kp (150) 

kp is the thermal conductivity of the porous media. The implicit 
method of finite-difference scheme had been employed to solve 
these models. The effect of the thermal conductivity of the porous 
media on the heat transfer characteristics and performance have 
been considered. The theoretical results obtained from the model 
were validated by comparison with the experimental data of 


Sopain et al. [70] and there was reasonable agreement. 
2.10. Double-pass flat plate solar air heater with longitudinal fins 


Naphon [71] has investigated heat transfer characteristics and 
entropy generation of the double-pass flat plate solar air heater 
with longitudinal fins. The basic physical equations used to 
describe the heat transfer characteristics of the solar air heater, as 
shown in Fig. 12a and b, are developed from the conservation 
equations of energy. 

With the first law of thermodynamics and based on the Naphon 
and Kongtragood’s model [17], Naphon developed a mathematical 
model for double-pass flat solar air heater with longitudinal fins. 
Using the assumption proposed in the case of double pass solar air 
heater with porous media, the following equations have been 
developed: 


For the glass cover: gS = hw(Tg — Ta) + he picg(Tg — T f1) 


+ Arg p(Tg — Tp) + hrag(Tg — Ta) (151) 
For the first-pass air : mC f- = he pig(Tg — T f1) 
N z=H 
+hefip(Tp-Tp) + f 2h (Ty — T pr)dz (152) 
frontal /z=0 
(a) Absorber i i 
Fin R hs) -l Glass cover 
hng Ph Air inlet 


Air outlet 


m) 
l a : i Absorber / 
hrp 


hz 


Fig. 12. (a) Schematic diagram of the solar heater (Naphon [71] model). (b) 
Schematic diagram of fin section (Naphon [71] model). 
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Fig. 13. (a) Three pass single cover air heater (Chandra et al. [72]). (b) Three pass two covers air heater (Choudhury et al. [75]). 


where m is the mass flow rate of fluid per unit with, Agontai is the 
frontal area of the channel, H is the height of the fin, z is the 
distance from the absorber plate, T, is the fin temperature in the 
upper channel and L is the length of the collector. N is the number 
of fins. 


N Z=H 
For the upper fin : i 2Lhy (Ty — T fı)dz 
frontal /z=0 
N dT 
= kA n) 153 
Afrontal ( dz z=0 ( ) 


For absorber plate: &gTgS = he fy p(T p — T f1) + hef2p(Tp — T f2) 


N dTy1 
Prep (Tp Ts) , Mrpo(Tp Th) 4 Afrontal ( aie dz Ja 
(a52) (154) 


T,2 is the fin temperature in lower channel and A, is the cross- 
section area of fin. 


For lower fin : 
( kA 


For second-pass air : 


z=H 
L 2Lhz(T,2 — T p2)dz 


z=0 


m 


frontal 


dT yp 
dz 


N 
Agrontal 


(155) 


dT p 
dx 
= hef2p(Tp — T f2) + he pop(T p -T f2) 


Z=H 
L 


For bottom plate : 
+ Ua(Tp — Ta) = 0 


2Lh2 (Tv2 = T f2)dz (156) 


Afrontal 


Ne f2b(Tpb — T f2) + hrpb(Tp — Tp) 
(157) 


It was found that the thermal efficiency increases with 
increasing the height and number of fins. 


2.11. Triple pass solar air heaters 


The schematic configuration of a typical triple-pass solar air 
heater is shown in Fig. 13a. It consists of one cover, an absorber 
plate and two rear plates (Fig. 13a). Chandra et al. [72] have 
investigated this design. Neglecting the thermal capacity of the 
plates, the energy balance equations were written as follows: 


Cover: gS = hrgs(Tg 
+ h;gpı (Tg = T f3) 


Ts) + hw(Tg — Ta) + heg f3(Tg — T f3) 


(158) 


Fluid3: -y E = lerlTs- T ys) — epi sa(Tp3—T 1) (159) 
Plate 1: TgQ pS = hig pi (Tg = Tp) + hepi f3(T p1 = T f3) 


+ hrpir2(Tp1 — Tp2) + hepi f2(Tp1 — T f2) (160) 


. —mC ç dT 
Fluid 2 : W FP = he p(T p =T f2) — Nepa f2(T f2 — T p2) 
161) 
Plate 2 : h;pıp2(Tp1 = T p2) T Ne f2p2(T f2 T Tp2) 
= hypap(T p2 — Ty) + hep2 pi (T p2 — T f1) 162) 
: mC dT 
Fluid 1 : w ra = hepfi(Tp2 —Tp)—ho p(T —Te) (163) 
Plate 3: hy pap (T p2 — T 3) + Rep p(T f — Tp) = U,(Tp = Ta) 164) 


with Ta(x = 0) = Ta, Tr (x = L) = Tja(L) and Tj2(0) = Ty3(0). 
In order to appreciate the analytical results, numerical 
calculations have been carried out corresponding to the experi- 
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mental data of Ezeike [74]. The results of this theoretical work 
differ substantially from Ezeike’s results. The results obtained from 
an exact analysis and from an approximate equivalent thermal 
network method were almost consistent. 

This model of solar air heater has also been investigated by 
Choudhury et al. [75]. They assumed that the air inlet temperature 
is the same as the ambient temperature. In the mathematical 
model proposed by these authors, Eq. (158) is reconsidered 
without the term of radiative exchange between the cover and the 
sky; Eqs. (160), (162) and (163) are the same. Eqs. (159), (161) and 
(163) have been replaced by the following: 


Fluid 3 : mC p(T fo = T f20) 
: = Neg 53(Tg — T f3) — hepi f3(T f3 — Tp) (165) 
Fluid 2: mC f(T 20 = T fio) 
o = hep f2(Tp — T f2) — Repo f2(T f2 — T p2) (166) 
Fluid 1: mC;(T fio — Ta 
= hep fı (Tp2 — T f1) — hep pa (T fi — Tp) (167) 


Another three-pass solar air heater was investigated by 
Choudhury et al. [75]. It consists of two covers, an absorber and 
a rear plate (Fig. 13b). Assuming the air inlet temperature to be the 
same as the ambient temperature, the steady state energy balance 
equations for the different plates and flowing air in the inlet, 
intermediate and outlet channels of the air heaters can be written 
as follows: 


Cover 1: AgS = hgia(Tgı = Ta) + h;g2g1 (Ta = Tg2) 
+ Nees fi (Tg1 — T F1) (168) 
Fluid 1 : mC ¢(T fio = Ta) 
= hegi fi (Tg1 = Tf) = heg2 f1 (Tg2 = Tf) (169) 
Cover 2 : Tg1Q p2S = Rega f1 (Tg2 = Tfi) T heg2 f2(Tg2 = T f2) 
+ hrg2g1(Tg2 = Tg) T hy pt p2(T pi = T p2) (170) 
Fluid2: CCT poo — T to) 
= Nega f2(Tg2 — T f2) — hepi f2(Tp1 — T f2) (171) 


(a) 


Absorber : Tg1 Tg2Q pS = hepi f2(T p1 = T f2) + Rept p(T pt = T f3) 
+ Ay pig2(T pr = Tg2) + Arpt p2(T pi = T p2) 172 
Fluid 3 : mC p(T fo = Ta) 
= hepi f3(T p1 —T 3) — Nepa f3(T f3 —Tp2) 173 
Rear plate: hy p1p2(T p1 — Tp2) + hcp2 f3(T f3 — T p2) 
= U;(T p2 — Ta) 174 


In all these equations Ta, Te and Tg are defined as 
Ta = (Trio + Ta)/2, Tp = (Treo + Thio)/2 and Tp = (Tf + Tpo)/2, respec- 
tively. In this study, the authors have made an effort to find the 
effects of duct length, duct depths of inlet (D1), intermediate (D2) 
and outlet (D3) air channel and mass flow rate (m) on the cost 
benefit ratio. Efforts were also made to identify the most superior 
air heater by comparing the cost effective performance of these air 
heaters with those of single pass with no, single and double cover 
and two pass with single and double covers air heaters studied 
earlier by the authors [33,75]. Choudhury et al. [75] in their 
analysis have used the heat transfer coefficients predicted by the 
equations given in Duffie and Beckmann [5], Tan and Charter [6] 
and Hollands and Shewen [61]. It was shown that the double cover, 
three pass solar air heater with air flow from top to bottom is 
observed to be more cost-effective than the single cover, three pass 
solar air heater with air flow from bottom to top. 


2.12. Multi-pass flat-plate solar air heaters with external recycle 


Ho et al. [76] have developed a theoretical formulation for a 
multi-pass solar air heater with external recycle and investigate 
the recycle effect on collector efficiency. The effects of channel 
thickness ratio on the collector efficiency was also studied. In this 
study, Ho et al. [76] have considered a recycled four-pass solar air 
heater with vertical thickness 2œH and 2(1 — @)H (a is the channel 
thickness ratio in vertical direction), and horizontal width 26H and 
2(1 — 6)H ($ is the channel thickness ratio in horizontal direction), 
respectively, which was obtained by inserting an absorbing plate 
and insulation sheet with negligible thickness into a parallel 


I7 
KEZ 


2(1-á4)H 


Ta(z) .Rm 


Absorbing, 


Fig. 14. (a) Multi-pass flat-plate solar air heater with recycle (Ho et al. [76]). (b) Schematic view of energy-flow of multi-pass type solar air heater with recycle (see Ho et al. 


[76]). 
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conduit of height 2H, length L, and width W (see Fig. 14). Assuming 
that: (i) the temperatures of the absorbing plate, bottom plate and 
bulk fluid were functions of the flow direction only, (ii) both glass 
covers and fluids do not absorb radiant energy, (iii) the 
absorptance of the glass cover 2 was negligible, the energy 
equations have been written as 


Cover1: Mrpei(T p — Tg1) + hy (Tc(X) — Tg1) + h4(Ta(x) — Tg1) 
gis(Tg1 D Ts) (175) 
Absorbing plate: &pTgıTg2S — hr(Tp — Ts) — hg(Tp — Ts) 
hy (Tp — Ta(x)) — ha(Tp — Ty (x)) + h3(Te(x) — Tp) + ha(Ta(x) 
_T,)=0 (176) 
Bottom plate: hy (Ta(x) — Tr) + h, (Tp (x) — Tr) + hrpr(Tp — Tr) 
= Ups(T p ~ Ts) (177) 
Channel a(flowing under the absorber plate) : 
hy (Tp — Ta(x)) — hy (Ta(x) — Tr) 
(R+1)mCp] dTq(x) 
=| WB dx (178) 
Channelb: hz(Tp-— Ty(x)) — hy (Tp (x) — Tr) 
(R+ 1)mC p] dT,(x) 
- [At BW Jx (179) 
Channel c: h3(Tp — Te(x)) — h3 (Te(xX) — Tg1) 
_ [-RmCp] dT (x) 
- Fa ao 
Channel d: h4(Tp — Ta(x)) — h4(Ta(x) — Tg1) 
_ | RmCp | dT)(x) 
awa i 
In this work, the equivalent diameters were defined as 
D. = (4a@H)W B [4aH|W(1 — £) 
ea (20H) + WB’ eb — [2aH] + W(1 — A)’ 
p.. - 4d = Hwee p, ,— 40 -o)HIWA — B) 
e 2A- aH] + WB’ ed (201 — aH] + WC — B) 
(182) 
and average air velocities given by the following: 
5, = Ma Mp 
~ QaH)WBp’ 5 2HaW(1 — pp’ 
Ue = ie J Ma (183) 
““ DHA —a)WBp’ “4 2H0 -a)W(1- bjo 


The above set of equations was solved numerically using an 
iterative procedure. In this procedure, the prescribed collector 
geometries, system properties and operating conditions were 
given, as well the initial guesses of the mean temperatures were 
specified for calculating the physical properties. The new values of 
mean temperatures were then re-calculated. If the re-calculated 
average temperature values were not quite close to the assumed 
values, the iteration calculations performed above will be done 
until the last assumed values meet the final calculated values. The 
results have shown that, the introduction of recycle effect was a 
feasible way to increase the collector coefficient due to the 
increasing fluid velocity effect and may compensate for the 
temperature difference decrement. 


2.13. Solar air heater with a compound parabolic concentrator 


Fig. 15 shows this solar air model. Pramuang and Exell [77] 
reported the results of an experimental study in which the method 
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Fig. 15. Schematic view of a solar air heater with a compound parabolic 
concentrator (Premuang and Exell [77] and Tchinda [78]). 


of Chungpaibulpatana and Exell is used to determine the collector 
parameters of a solar flat plate collector with a CPC for heating air. 
Recently, Tchinda [78] has quantified the heat transfer within this 
design of solar air heater and a mathematical model analysing the 
collector thermal performance has been proposed. In this model 
(see Fig. 15), the following assumptions were made: (i) It was 
assumed that the CPC is ideal and free from fabrication errors. (ii) 
Any beam of radiation incident within the acceptance angle 6,, 
with the help of the parabolic reflector can reach the receiver. (iii) 
The reflection of radiation from the parabolic reflector was 
accounted for by the apparent reflectance p,<"* . (iv) The direction 
of the beam radiation incident on various components in the 
collector can be found though geometry. (v) The succeeding 
absorption and transmission processes inside CPC were diffusive 
and were accounted for in terms of diffuse properties. The solar and 
infrared energy exchanges in the collector were treated separately 
using pertinent radiative properties in the spectrum. (vi) Physical 
and optical properties of materials were assumed to be indepen- 
dent of temperature. (vii) The concentrator did not produce an 
image of the light source; hence it was called non-imaging 
concentrators. With the above assumptions, the energy equations 
may be written as follows: 


T ee A 
For cover : MM pe SE =I(t)[ae+ GTP phn” | a + hrp(Tp — Te) 
p 
t hp/c(Tp Tc) — hrs(Te — Ts) Neja(Tec Tp) (184) 
with t > 0. 
aT p 
For plate absorber: MpCpp 7a 
KOZ o Ap] Ac 
=I(t)TcPm' P| Ap 4 A pPpPc A” Ap hrp(Tp — Tc) Apjc(Tp Te) 
— qu (t) (185) 
with t > 0. 
For fluid : pei! 
MC yf OT 
dult) — PE aye = Uo(T y ~ To) (186) 


with t>0 and 0 <x < L. where P=1 — g/l, ([80]), g is the gap 
thickness. A. = W*L, Ap=I,"L and qu„= UT,- T). An iterative 
method was used by the author to bring the effect of the 
temperature dependence of various heat transfer coefficients. For 
certain temperatures, they were first calculated by using the 
standard expressions given earlier. Equations were solved by 
assuming the constant and the new solutions used to generate all 
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Fig. 16. (a) v-Groove air collector, Karim and Hawlader [83]. (b) v-Groove air collector 
v-Groove air collector (Kabeel and Mecarik [50]). 


the heat transfer coefficients again till the values converge. This 
model has been validated with the experimental data of Pramuang 
and Exell [77]. 


2.14. V-groove solar air collector 


The design given in Fig. 16a was also studied experimentally 
and theoretically by Karim and Hawlader [82,83]. The energy 
balances for absorber plate, back insulation and working fluid of 
the above collector are given by the following Karim and Hawlader 
[83]: 


TgdpS + U:(Ta — Tp) +h (Tp — Tp) thy Tp) = 0 (187) 
Up(Ta — Ty) + h2 (T f — Ty) + hr(Ty — Tp) = 0 (188) 
hy(Tp —T >) +ho(Tp — Tp) = qu (189) 


A series of experiments were conducted, data were recorded for 
different operating variables to determine the performance of v- 
groove solar air collectors. Experimental and analytical results 
have shown a good thermal performance of the v-groove collector. 
Satisfactory qualitative and quantitative agreement between 
experimental and analytical results was achieved. 

An Internet based mathematical simulation has been conducted 
and developed by Bashria et al. [84,85] and Bashria et al. [86] to 
predict the thermal performance for different designs of solar air 
heaters with v-groove absorbers, single, double or triple glass 
covers. 

Fig. 16a-c show the designs of v-groove solar air collector 
studied. A number of simplifying assumptions has been made to 
lay the foundation without obscuring the basic physical 
situation. These assumptions are as follows: (i) performance is 
steady state; (ii) There is no absorption of solar energy by a cover 
insofar as affects losses from the collector, (iii) Heat transfer fluid 
is considered a non-participating medium, (iv) The radiation 


(d) 


S 
| upper glass cover 


Back insulation 


Up 


Bashria et al. [84]). (c) v-Groove air collector with porous media (Bashria et al. [85]). (d) 


coefficient between the two air duct surfaces is found by 
assuming a mean radiant temperature equal to the mean fluid 
temperature. (v) Loss through front and back are of the same 
temperature. The steady state energy equations yield the 
following equations: 

In the single pass v-groove absorber: 


For the cover: U;(Tg — Ta) = hı (Tp — Tg) + hr (Tp — To) (190) 


v-groove absorber: Tg@pS = hy (Tp — Tg) + hr(Tp — Tg) 
+ ho(Tp — Tf) (191) 


CA dT 
Fluidmedium: h(Tp — T f) + Up(Ta T- (A) Ze (192) 


In double pass double duct: 


Forthe cover: U;(Tg — Ta) = hi (Tp — Tg) + hy (Tp — Tg) 193 


v-groove absorber: Tg pS = h2(Tp — T 1) + Wi (Tp — Tg) 


t hy(Tp — Tg) + h3(Tp — T f2) + h,2(Tp — Tr) 194 
Fluid medium in the upper passage: —h2(T p — T 1) 
_ mC p dT f 
n (o ) z 195 


Fluid medium in the lower duct: h3(Tp — T f2) + h4(T f2 — Tr) 
= (WW) dx une 
Bottom plate: Up(T; — Ta) = ha(T f2 — Tr) + h2 (Tp — Tr) (197) 


In double flow double duct packed: 


For the cover: U;(Tg — Ta) = hı (Tp — Tg) + hr(Tp — Tz) (198) 
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v-groove absorber: TgQ@pS = h2(Tp — T p,) + hi (Tp — Tg) 


+ hy(Tp — Tg) + h3(Tp — T f2) + h2(Tp — T pr) (199) 
Fluid medium in the upper passage: (Tp — T f1) 
B MC p dT fi 
- (me) 200 


Fluid medium in the lower duct : 


mC dT 
+T- T= (TO G2 


h3(Tp — T f2) + h4(T f2 — T pr) 
(201) 


Porous media: h4(Tpr — T f2) = h2(Tp — Tpr) + he(T pr — Tr) 
(202) 
Bottom plate: Us(T f2 — Tr) + he(T pr — Tr) = Up(Tr — Ta) (203) 


The pressure drop through the channel in the air heater has 
been computed using the following expression [87]: 


D] /m?\ /L\?3 
P= [ntr] (%) 0) 
where for fo=24/Re, y=0.9 for laminar flow (Re < 2550), 
fo = 0.0094, y=2.92Re°!° for transitional flow (2550 < Re < 
10*) and f,=0.059Re°*, y=0.73 for turbulent flow 
(104 < Re < 10°). 

This mathematical model has been used recently by Bashria 
et al. [86] to find the effect of different parameters such as mass 
flow rate, flow channel depth and collector length on the system 
thermal performance and pressure drop through the collector (see 
Fig. 16a—c). The prototype internet program was developed to be 
used as a tool to support the design of solar air heaters used for 
drying. The results were validated by comparison the predicted 
output and the experimental results carried out by Karim and 
Hawlader [82]. A great correlation has occurred between the 
experimental and the predicted efficiencies. Also a good correla- 
tion has been found between the experimental and predicted 
efficiencies and outlet temperature, respectively. 

Kabeel and Mecarik [50] have investigated v-groove air heater 
system with two covers and the flow passes over and under the 
absorber plate (Fig. 16d). By considering the energy balance for 
each part of the system in the steady state heat flow were written 
as the follow: 


(204) 


For the upper glass cover: @&@g1S + Nygigo(Tg2 — Tg1) 
+ hegia (Tar — T gi) — hw(Tg1 — Ta) = 0 205) 


For the second glass cover: Tg1Q@g2S — heg2a1(Tg2 — Tar) 
= ho(T f2 — Tg2) — hrg2g1 (Tg2 — Tg1) + hrg2p(Tp — Tg2)=0 (206) 


For the flow between the two glass covers: hg2a1(Tg2 — Ta1) 
= hga (Ta = Tg1) =0 207) 


For the absorber plate : 
— Zhpf2(Tp — T f2) — hrpg2(Tp — Tg2) 


Tg1 Tg2 pS — ZNp fi (T p — T fi) 
hrpo(Tp —Tp) =0 (208) 


For Back plate: hy pp(Tp — Ty) + Reb pi (T fı — Th) — Up(Tp — Ta) 
=0 209) 


Z is the ratio between the actual absorber area and the absorber 
area normal to the solar radiation following upon the collector. 
With this mathematical model taking into account the absorber 
shape factor, Kabeel and Mecarik [50] have outlined the effect of 
the change of the triangular angle 6 on the pressure drop, on the 
heat transfer added to the flowing air and deduced the optimum 
angle of the v-groove. 


2.15. Single pass solar air heater with packed 


2.15.1. Type 1 

The type 1 solar heater is shown in Fig. 17a. It consists of 
conventional configuration of Fig. 4 with the flow passage filled 
with packing. In the mathematical model proposed by Choudhury 
and Garg [88], energy balance equation of the cover is given by 
Eq. (36). For the absorber plate, packing and air flow and bottom 
plate, energy balance equations in the steady state were given by 
the follow: 


Absorber plate: Tg@pS = (hrgp + Ncgp)(T p — Tg) 


+ hrppi (Tp — Tp) + Hep s(Tp — Tf) (210) 


Packing: = Mrppi(Tp — Tpi) = hepi p(T p1 -T f) + hrpib(Tpi — To) 
(211) 
: mC , dT 
Air flow: 4 EF = heps(Tp — Tp) + hep p(T pt — Ty) 
+ Ney p(T — Tf) (212) 


Bottom plate: h,pıb(Tpi — Th) = Aen f(T, — T f) + Up (Tp — Ta) 


(213) 


The boundary conditions for the above design are: 
T(x = 0) =Tg, T(x =L) =T. In this model, it was assumed that 
the heat transfer resulting from conduction in the packed flow 
channels is neglected due to point contacts between packing and 
the plates. 

Recently, Paul and Saini [90] have used the mathematical model 
proposed by Choudhury and Garg [88] to outline the bed 
parameters in such a way that the solar energy collection systems 
deliver energy with minimum cost. A computer program in ‘C 
language has been developed to obtain thermal efficiency and cost 
per unit energy delivered by packed bed solar air heater having 
wire mesh screens and pebbles as packing materials with different 
specifications. This program was validated by comparing some of 
the predicted values of thermal efficiency with experimental 
values obtained by Collier [92]. 


2.15.2. Type 2 

It consists of two sheets of cover and the back plate; the 
passage between the inner cover and the back plate is provided 
with black-painted packing, which act as absorbers of solar 
radiation and transfer the heat to the air flowing through the 
packed channel (see Fig. 17b). This device has been investigated 
by Choudhury and Garg [88] and the energy balance written. 
Energy balance of the outlet cover, back plate and air flow 
were given by Eqs. (36), (212) and (213) respectively. For 
the inner cover and packing, the energy balances are as 
follows: 


For the inner cover: Tg2Œg1S + hy pig2(T pı — T p2) 


= (h,g1g2 + hegig2)(Tg2 = Tg ) T heg2f(Tg2 = Tj) (214) 
For the packing: TeffÆpS = hepi f(T pi + Tp) 
= hrpig2(Tp1 + Tp2) 
+ h,pib(T pi + Tp) (215) 


This design has been used as a component of hybrid heating 
system [93,94]. To simulate the performance of this system, a 
transient analysis has been proposed by the authors where the 
above energy balance has been written for a rock bed solar air 
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Fig. 17. (a) A schematic view of a single cover solar air heater (Choudhury and Garg [88]). (b) A schematic view of a single cover solar air heater (Choudhury et al. [88]). 


heater. 
dT 
Mei Cg at ~ AgS + Ngig2(Tg2 —T 1) — Agta(Te1 — Ta) (216) 
dis 
Mg2Ce2 -e = TgiQ@g2S + Mpiga(T p1 — Tg2) — Ngagi (Tg2 — Tg1) 
— hg2a(Tg2 -T f) (217) 
dT 
MC a = Tg1Tg2 pS = hepi f(T p1 — Tp) 
— h; pıg2(T p1 — Tg2) — hpib(T pi — To) (218) 
dT ; 
M pC p—ae = CT- Tp) + Apr f(T pr — T 5) — Nga p(Te2 
— Ty) — hey (Tp — T 5) (219) 
dT, 
MoCo -iF = hpb(Tp1i — Tb) — Neb p(T — Tf) — Ub(Tp — Ta) (220) 


The results have shown that the hybrid air-to-water heating 
system performs much better when coupled with a rock bed air 
heater than when coupled to a conventional empty channel air- 
heating collector. 


2.16. Double-pass solar air heater with packed 


The schematic diagrams of the double-pass solar air heater with 
packed bed are shown in Fig. 18. The different heat transfer 
mechanisms in terms of the various heat transfer coefficients are 
also shown. The air was firstly forced through the packed bed 
existing in the upper channel (Fig. 18a), formed between the lower 
cover and the absorber plate, and is then re-circulated to flow in 
the opposite direction through the lower channel, formed between 
the absorber and back plates. Ramadan et al. [96] have investigated 
heat transfer in this solar heater. The mathematical model has been 
developed taking account of the following assumptions: (i) the air 
heater operates under steady state conditions. However, the model 


remains transient due to the time dependence of solar radiation 
intensity and ambient temperature. (ii) The heat capacities of the 
glass covers, absorber and back plates and insulation were 
negligible. (iii) The temperature of the flowing air was varied 
only in the direction of flow. (iv) The thermo physical properties of 
the flowing air were assumed to vary linearly with temperature 
(see [7]). (v) There was no temperature gradient across the 
thickness of the covers, the absorber and the back plates. On the 
basis of the above assumptions, the energy balance equations for 
the different components of the air heater were written as follows: 


Cover2: Ag@gS = Ag(Hrgig2 + Negig2)(Tg2 — Tg1) 
+ hwAg(Tg2 = Ta) Eg Aghrgasky(Tg2 = Tsky) 
(221) 
Cover 1: AgQ@gTgS = Ag(ħrgig2 + Ncgig2)(Te1 — Tg2) 
= h; pig1Ag(T p1 = Ta) = Agheei f1 (Tf = Tg) (222) 


Air flow in upper channel: Whpfi(Tp-T f) 


dT 
+Whpfi(Tp -Tf)= macs + Whey f1(T fi — Tgi) 


+WU,(T f — Ta) 223 


Absorber plate: Aphrp1p(Tp1 — Tp) = Aphp pi (Tp — Tfi) 
+ ApNp fa(T p — T f2) + Aphrpo(T p — To) 224 


Air flow in lower channel: Why f2(T p — T f2) 


T 
P $ WUG(T p — Ta) 225 


d 
+ Wh f(T -T pr) = My2C fo dx 


Back plate : Abhy f2 (T f2 — Tp) + Abh; pb(Tp — Tp) 
= ApUp(Tp — Ta) 226 


The various heat transfer coefficients given in the above 
equations were calculated using the correlation given in the 
literature [7]. A suitable computer program was prepared, for the 
solution of the energy balance equations for different elements of 
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Fig. 18. (a) Schematic view of energy-flow of double-pass type solar air heater with recycle (see Ramadan et al. [96]). (b) Schematic view of energy-flow of double-pass type 


solar air heater with recycle (see El-Sebaii et al. [99]). 


the solar air heater. It was inferred that for increasing the outlet 
temperature of the flowing air after sunset, it was advisable to use 
the packed bed materials with higher masses and therefore with 
low porosities. It was recommended to operate the system with 
packed bed with values of the mass flow rate equal 0.05 kg/s or 
lower to have a lower pressure drop across the system. This 
mathematical model has been validated by comparisons between 
experimental and theoretical results showing that good agreement 
was achieved. 

El-Sebaii et al. [99] have proposed an investigation of heat 
transfer on a double-pass solar air heater in which the air was 
firstly forced through the upper channel, formed between the 
lower cover and the absorber plate, and then re-circulated to 
flow in the opposite direction through the packed bed that exists 
in the lower channel, formed between the absorber and back 
plates (Fig. 18b). With the above assumptions, El-Sebaii et al. 
[99] have written the energy balance equations for the system. 
For cover 2, energy balance equation is given by Eq. (221). For 
cover 1, air flowing in upper channel, absorber plate, air flowing 
in the lower channel and back plate may be written as follows, 
respectively: 


AgQtgT gS — Ag(Mrgiga + hcgig2)(Tg1 — Tg2) 


= —hrpigiAg(T pt — Tg1) — Aghtegi pi (T fı — Tg1) (227) 
Air flow in upper channel: Why i (Tp — Tp) 
dT 
+ Whey pi (Tp — Tp) = mpc + Wheg1 f1(T pi — Tgi) 
+ WUS(T fi — Ta) (228) 


Absorber plate: = @ptgApS = Aphep p(T p — T f1) 

+ Aphep fa(T p — T f2) + Aphrpgi (T p — Tg1) + Aphrppi(T p — Tp) 

(229) 
Air flow in lower channel: Why f2(T p -T f2) 
dT 

+ Whepr sa(T pt — Tyo) = M pa f2 = + WUS(T pa — Ta) 

+ Wh; p(T f2 — Tp) (230) 
Back plate : Aphp f2(T f2 — Tp) + Aphy pip (T p1 — Tp) 

= ApU, (Tp — Ta) (231) 


The various heat transfer coefficients given in the above 
equations were calculated using the correlation given in literature 
[97,98]. The program developed by Ramadan et al. [96] was used. 
On the basis of the obtained results for the solar air heater, it had 
been concluded that (i) to enhance the convective heat transfer 
coefficients to the flowing air, as well as the heater efficiency, it 
was advisable to operate the system with the packed bed above the 
absorber plate with low porosity and relatively high mass flow 
rates of air. (ii) It was better to use gravel instead of limestone as a 
packed bed above the absorber plate. (iii) It was recommended to 
operate the system with and without the packed bed with values of 
the mass flow rate equal to 0.05 kg/s or lower to have a lower 
pressure drop across the system, and therefore, reasonably high 
thermo hydraulic efficiency. (iv) It was better to use materials with 
high thermal conductivity, such as iron scraps, as packed beds 
above the absorber plate if a higher outlet temperature of air is 
required during sunshine hours. Gravel as a packed bed above the 
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Fig. 19. Three pass two covers air heater with in-built thermal storage (Jain and Jain 
[100], Jain [101]). 


absorber plate during low intensity solar radiation periods was 
recommended. 


2.17. Multi-pass solar air heater with in-built thermal storage 


Recently, Jain and Jain [100] and Jain [101], studied an indirect 
solar crop dryer coupled with the solar air heater, investigated the 
heat transfer in the multi-pass solar air heater with in-build 
thermal storage. The schematic diagram is shown in Fig. 19. In this 
model, the solar radiation transmits from the glass covers and is 
absorbed by the absorber plate. The air flows in between the 
covers, above the absorber plate and below the storage material, 
where it is heated along the path. Jain and Jain [100] and Jain [101] 
have written the energy balance equations for this design on the 
base of the following assumptions (i) the heat capacities of the air, 
covers, absorber plate and insulation were negligible; (ii) there was 
no temperature gradient along the thickness of cover; (iii) storage 
material has an average temperature (T,) at a time t, (iv) there was 
no stratification existing perpendicular to the air flow in ducts, (v) 
the system was perfectly insulated and there was no air leakage. 
Based on these assumptions, energy balance equations on the 
double-cover multi-pass flat plate collector with in-built thermal 
storage were given as follows: 


For the first cover: &Œg1Ag1ıSeff + Nrgagi1Agi(Tg2 — Tg1) 
= Negi ptAgi (Ta = Tf) ae AegiaAgi (Tg1 z Ta) 


ste h;giskyAg1 (Tg1 = Tsky) (232) 
where Tsky = Ta — 6 (see [102]): 
For the second cover: gz Tg1Ag2Se f 5 + Mrpg2Aga(T p 

a Tg2) + he f1g24g2(T fı = Tg2) = Aygag1 Aga (Tg2 = Tg1) 

+ hogz f2Ag2(Tg2 — T f2) (233) 


For the absorber plate: Œ pTg1ıTg2Ag2Seff = hrpg2Ap(Tp — Tg2) 


+ hepf2Ap(Tp — T f2) + hcpsAp(Tp — Ts) (234) 
For air flow in the first channel: Wheg 51 (Tg1 — T f1) 

= mC E + Whe figa(T fi — Tg2) (235) 
For air flow in the second channel: Wheg2f2(Tg2 — T f2) 

=mCy = + Whepj2(T f1 — Tg2) (236) 
For air flow in the third channel: Whea p3(Ts — T f3) 

=mCy a + WU4(T 3 — Ta) (237) 
For the storage material: hcps(Tp — Ts)Ap 

= M,C; ae + Nes p3(Ts — T p3)Ap (238) 


The convective heat transfer coefficients from the plate to cover, 
parallel to each other and inclined at the angle £ to horizontal has 
been expressed given in literature [7,27] and in the appendix. Jain 


and Jain [100] have computed solar intensity on the inclined 
surface by using the method given by Lui and Jordan [103]. A 
computer program was prepared to solve the above set of 
equations. The performance evaluation of a tilted multi-pass solar 
air heater with in-built thermal storage has been carried out for 
deep-bed drying applications. The thermal energy storage also 
affect during the off-sunshine hours is very pertinent for crop 
drying applications. With the above mathematical model, Jain 
[101] has shown that it is useful to evaluate the thermal 
performance of a flat plate solar air heater for the crop drying in 
multiple trays. 


2.18. Matrix (porous) solar air heater 


The porous or matrix type solar air heaters are shown in Fig. 20a 
and b. It consists of single glazing at the top, matrix, back metallic 
plate and insulation. In the first design (Fig. 20a) the air flow is 
upward through the matrix and in the second design (Fig. 20b), it is 
flow downward through the matrix. These models have been 
investigated by Sharma et al. [98] and Sharma et al. [104]. Sharma 
et al. [98] analyzed the various energy exchanges inside the porous 
bed. The energy balance equations were described by the following 
equations for air flowing upward: 


aT FT 

M,C, mi kp ay? h p(T, — T f) + waS exp(—ay) (239) 
aT OTs 

MCs- = GCy By h y(Tp — Tf) (240) 


The corresponding relations for air flowing downward can be 
obtained from the above equation by changing the sign of G. The 
above equations were solved explicitly. The results were showed 
that such types of solar air heater have the tendency to maintain a 
constant efficiency throughout the operation. The experiment was 
also carried out. A good agreement of the theory with experimental 
observations almost confirms that the model can be used to predict 
the performance of the porous solar air heaters. 

A transient model for a matrix air heater having a single glazing 
at the top and a thin metallic plate backed by an insulator at the 
bottom was presented by Sharma et al. [104]. The mathematical 
model used was similar to the one obtained by Sharma et al. [98] in 


“| Air outlet 


Back metallic 


insulation Plate 


Fig. 20. (a) Matrix air heater (air flow upward) (Sharma et al. [98]). (b) Matrix air 
heater (air flow downward) (Sharma et al. [98,104]). 
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the air flow downward air heater. In this model solar insulation and 
ambient temperature were represented by Fourier series in time 
with y dependent Fourier coefficients. The experimental results for 
a varying inlet air flow rate were analyzed and the rating 
parameters were included. From the results, it was observed that 
the matrix air heater concept shows promise for collection 
enhancement of energy at a cost lower than most of other 
conventional solar air heaters. 


3. Conclusion 


Many mathematical approaches have been developed for 
simulating the performance of solar air heaters. In this work, a 
review is given of the present state of mathematical modeling 
based on the energy analysis. 

The models on energy analysis are usually expressed in 
relatively simple equations and are based on the steady state of 


2(T pm/520) 


where the units of hy and V are W/m? C and m/s, respectively. If the 
unit of hw is changed to Kj/h m? K, the above equation becomes: 


hw =4.9+3.2V (A2) 


This coefficient has also been given by Watmuff, Charters and 
Proctor [109] as: 


hy =2.8+3.0V (0<V<7m/s) (A3) 


An empirical equation for Ur was developed by Klein [67] 
following the basic procedure of Hottel and Woertz [68]. For a 
horizontal collector with two glass covers Ur is given by the 
relationship: 


Ur = Ur + Ur2 
With 


(A4) 


U. pa 
ý [r + (1 +0.089)hw — 0.1166hwep 


={ 
(1 + 0.7866 x gera] 


o(Tpm + Tsky) (Tm + Tiy) 


Ur2 


the system. Certain assumptions used in the mathematical form 
appear to be too ideal. However these models are limited by the 
heat transfer efficiency between the absorber and the fluid. Since 
heat transfer augmentation usually increases friction, optimum 
parameters, such as mass flow and collector geometry, have to be 
determined. 

Model validation is a key step in model development since it 
offers the possibility of comparing computed results with actual 
system behavior. Experiments and/or comparison with known 
analytical, or previous numerical results are mostly used to 
validate the mathematical model. Good agreement is generally 
observed and the computed trends are also compatible with those 
reported for experimental observations. However, there appears to 
be a lack of agreement among authors in several areas and 
furthermore, there are still some key deficiencies that can only be 
addressed through continued enhancement of the existing 
modeling techniques. As such, there exists a large body of work 
to predict several features of the process. 

Future work should be directed towards a more comprehensive 
description of heat transfer phenomena, which would involve the 
generation of more information in the exergy analysis of the solar 
air heaters. 
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Appendix A 


The convective heat transfer coefficient họ for air flowing over 
the outside surface of the cover depends primarily on the wind 
velocity V. MacAdams [19] obtained the experimental result: 
hy = 5.7+3.8 V 


(0<V<5 m/s) (A1) 


~ (Ep +2 x 0.00591 hw) | + [2 x 2 + (1 + 0.089 hw — 0.1166 hw&p)(1 + 0.07866 x 2) — 1 + 0.133 €p] /£p — 2 


References 


= 


Chandra R, Sodha MS. Testing procedures for solar air heaters: a review. 
Energy Conversion and Management 1991;32(1):11-33. 

2] Shaobo H, Zeng D, Ye S, Zhang H. Exergy analysis of the solar multi-effect 
humidification-dehumidification desalination process. Desalination 2007; 
203:403-9. 

3] Ekechukwu OV, Norton B. Review of solar-energy drying systems III: low 
temperature air-heating solar collectors for crop drying applications. Energy 
Conversion and Management 1999;40:657-67. 

4] Choudhury C, Chauhan PM, Garg HP. Design curves for conventional solar air 
heaters. Renewable energy 1995;6(7):739-49. 

5] Duffie JA, Beckmann WA. Solar engineering of thermal processes. New York: 
Wiley; 1980. 

6] Tan HM, Charters WS. Effect of thermal entrance region on turbulent forced 
convective heat transfer for an asymmetrically heater rectangular duct with 
uniform heat flux. Solar Energy 1969;12:513. 

7] Ong KS. Thermal performance of solar air heaters: mathematical model and 
solution procedure. Solar Energy 1995;55(2):93-109. 

8] Njomo D. Unglazed selective absorber solar air collector: heat exchange 
analysis. Heat Mass Transfer Warme 2000;36:313-7. 

9] Njomo D, Daguenet M. Sensitivity analysis of thermal performances of flat 
plate solar air heaters. Heat Mass Transfer 2006;42:1065-81. 

10] Sparrow EM, Ramsey JW, Mass EA. Effect of finite width on heat transfer and 
fluid flow about an inclined rectangular plate. Journal of Heat Transfer— 
Transactions of the ASME 1979;101:199-204. 

11] Mercer WE, Pearce WM, Hichcock JE. Laminar forced convection in the 
entrance region between parallel flat plates. Journal of Heat Transfer— 
Transactions of the ASME 1967;89:251-7. 

12] Clark D. Passive/hybrid comfort cooling by thermal radiation. In: Passive/ 
hybrid cooling conference, Miami; 1981. 

13] Garg HP, Chandra R, Rani U. Transient analysis of solar air heaters using finite 
differences technique. Energy Research 1981;5:243-52. 

14] Njomo D. Modelling the heat exchanges in a solar air heater with a cover 
patially transparent to infrared radiation. Energy Conversion and Manage- 
ment 1991;31(5):495-503. 

15] Mohamad AA. High efficiency solar air heater. Solar Energy 1997;60(2): 
71-6. 

16] Hegazy AA. Performance of the flat solar air heaters with optimum channel 
geometry for constant/variable flow operation. Energy Conversion and Man- 
agement 2000;41:401-17. 

17] Paisarn N, Kongtragool B. Theoretical study on heat transfer characteristics 
and performance of the flat-plate solar air heaters. International Commu- 
nications in Heat and Mass Transfer 2003;30(8):1125-36. 

18] Zhai XQ, Dai YJ, Wang RZ. Comparison of heating and natural ventilation in a 
solar house induced by two roof solar collectors. Applied Thermal Engineer- 
ing 2005;25:741-57. 

19] McAdams WH. Heat transmission. New York: McGraw-Hill; 1954. 


1758 R. Tchinda/Renewable and Sustainable Energy Reviews 13 (2009) 1734-1759 


20] Kreith F, Kreider JF. Principles of solar engineering. New York: McGraw-Hill; 
1978. 

21] Heaton HS, Reynolds WC, Kays WM. Heat transfer in annular passages. 

Simultaneous development of velocity and temperature fields in laminar 

flow. International Journal of Heat and Mass Transfer 1964;7:763-81. 

22] Hausen H. Darstellung des warmenberganges in rohren durch verallgemei- 

nerte potenzbeiziehungen. VDIZ 1943;4:91-8. 

23] Tan HM, Charters WS. an experimental investigation of forced-convective 

heat transfer for fully-developed turbulent flow in a rectangular duct with 

asymmetric heating. Solar Energy 1970;13:121-5. 

24] Niles PW, Carnegie EJ, Pohl JG, Cherne JM. Design and performance of an air 

collector for industrial crop dehydration. Solar Energy 1979;20(1):19-23. 

25] Yeh HM, Lin TC. The effect of collector aspect ratio on the collector efficiency 

of flat-plate solar air heaters. Energy 1995;20(10):1041-7. 

26] Hegazy AA. Thermohydraulic performance of heating solar collectors with 

variable width, flat absorber plates. Energy Conversion and Management 

2000;41:1361-78. 

27] Aboul-Enein S, El-Sebaii AA, Ramadan MRI, El-Gohary HG. Parametric study 

of a solar air heater with and without thermal storage for solar drying 

applications. Renewable Energy 2000;21:505-22. 

28] Hollands KGT, Unny TE, Raithby GR, Konicek L. Free convective heat transfer 

across inclined air layers. Journal of Heat Transfer—Transactions of the ASME 

1976;98(2):189-93. 

29] Siegel R, Howel JR. Thermal radiation heat transfer. Washington: Hemisphere 

Publishing Corporation; 1981. 

30] Njomo D. Techno-economic analysis of a plastic cover solar air heater. Energy 

Conversion and Management 1995;36(10):1023-9. 

31] Njomo D. Etude théorique du comportement thermique d’un capteur solaire 

plan a air a couverture combine plastique vitre. Revue Generale de Thermi- 

que 1998;37:973-80. 

32] Garg HP, Datta G, Bhargava K. Some studies on the flow passage dimension for 

solar air testing collector. Energy Conversion and Management 1984;24(3): 

181-4. 

33] Choudhury C, Chauhan PM, Garg HP. Performance and cost analysis of two- 

pass solar air heaters. Heat Recovery Systems & CHP 1995;15(8):755-73. 

34] Al-Kamil MT, Al-Ghareeb AA. Effect of thermal radiation inside solar air 

heaters. Energy Conversion and Management 1997;38(14):1451-8. 

35] Jannot Y, Coulibaly Y. Radiative heat transfer in a solar air heater covered with 

a plastic film. Solar Energy 1997;60(1):35-40. 

36] Hegazy AA. Optimum channel geometry for solar air heaters of conventional 

design and constant flow operation. Energy Conversion and Management 

1999;40:757-74. 

37] Hollands KGT, Shewen EC. Proc. Int. Sol. Energy Soc. Congr., SUN II, Silver 

Jubilee Congr, Atlanta, GA; 1979. 

38] Biehl FA. Test results and analysis of a convective loop solar air collector. In. 

Proceedings of 6th National Passive Solar Conference Portland, Oregon (cited 

by Alfeld K et al., 1988); 1981. 

39] Jha RK, Choudhury C, Garg HP, Zaidi ZH. Performance prediction of a solar 

heater house. Energy Conversion and Management 1992;33(4):263-73. 

40] Ucar A, Inalli M. Thermal and exergy analysis of solar air collectors with 

passive augmentation technique. International Communication in Heat and 

Mass Transfer 2006;33:1281-90. 

41] Wong HY, Handbook of essential formula and data on heat transfer for 

engineers, London; 1977. 

42] Siedel EN, Tate GE. Heat transfer and pressure drop of liquids in tubes. 

Industrial and Engineering Chemistry 1936;28:1429-35. 

43] Forson FK, Nazha MAA, et Rajakaruna H. Experimental and simulation studies 

ona single pass, double duct solar air heater. Energy Conversion and Manage- 

ment 2003;44:1209-27. 

44] Achenbach PR, Cole SD. Performance of fourteen masonry chimneys under 

steady state conditions. Transactions of the ASHVE 1949;55:129-54. 

45] Yeh Ho-Ming, Chii-Dong Ho, Jun-Ze Hou. The improvement of collector 

efficiency in solar air heaters by simultaneously air flow over and under 

the absorbing plate. Energy 1999;24:857-71. 

46] Kays WM. Convective heat and mass transfer. New York: McGraw-Hill; 1980. 

47] Yeh Ho-Ming, Ho Chii-Dong, Hou Jun-Ze. Collector efficiency of double-flow 

solar air heaters with fin attached. Energy 2002;27:715-27. 

48] Al-Nimr MA, Damseh RA. Dynamic behaviour of baffled solar air heaters. 

Renewable Energy 1998;13(2):153-63. 

49] Matrawy KK. Theoretical analysis for an air heater with a box-type absorber. 

Solar Energy 1998;63(3):191-8. 

50] Kabeel AE, Mecarik KK. Shape optimization for absorber plates of solar air 

collector. Renewable Energy 1998;13(1):121-31. 

51] Ammari HD. A mathematical model of thermal performance of a solar air 

heater with slats. Renewable Energy 2003;28:1597-615. 

52] Duffie JA, Beckmann WA. Solar engineering of thermal processes, 2nd ed., 

Toronto: John Wiley & Sons, Inc.; 1991. p. 949. 

53] Nusselt W. Der warmeaustatsch zwischen wand und wasser im rohr. For- 

shung auf dem Gebiete des Ingenieurwesens 1931;2:309. 

54] Swinbank WC. Long-wave radiation from clear skies. QJR Meteorological 

Society 1963;89-339. 

55] Satcunanathan, Deonarine. A two-pass solar air heater. Solar Energy 

1973;15:41. 

56] Caouris Y, et al. A novel solar collector. Solar Energy 1978;21:157. 

57] Wijeysundera NE, Lee Lee AH, Tjioe LK. Thermal performance study of two- 
pass solar air heaters. Solar Energy 1982;28(5):363-70. 


[58] Wijeysundera NE. A net radiation method for the transmittance and absorp- 
tivity of series of parallel regions. Solar Energy 1975;17:75. 

[59] Verma Ratna, Chandra Ram, Garg HP. Optimization of solar air heaters of 
different designs. Renewable Energy 1992;2(4/5):521-31. 

[60] Choudhury S C., Garg HP. Comparative performance of Coriander dryer 
coupled to solar air heater and solar air-heater-cum-rock bed storage. 
Applied Thermal Engineering 1996;16(6):475-86. 

[61] Hollands KGT, Shewen EC. Optimization of flow passage geometry for air- 
handling plate type solar collectors. Journal of Solar Energy Engineering 
1981;103:323-30. 

[62] Hirunlabh J, Kongduang W, Namprakai P, Khedari J. study of natural ventila- 
tion of houses by a solar wall under tropical climate. Renewable Energy 
1999;18:109-19. 

[63] Yamali Cemil, Ismail Solmus. Theoretical investigation of humidification- 
deshumidification desalination system configured by a double-pass flat solar 
air heater. Desalination 2007;205:163-77. 

[64] Incropera FP, Dewitt DP. Fundamentals of heat and mass transfer, 3rd ed., 
Wiley; 1990. 

[65] Yeh Ho-Ming, Chii-Dong Ho, Wen-Song Sheu. Double-pass heat or mass 
transfer through a parallel-plate channel with recycle. International Journal 
of Heat and Mass Transfer 2000;43:487-91. 

[66] Ho CD, Yeh HM, Wang RC. Heat-transfer enhancement in double-pass flat- 
plate solar air heaters with recycle. Energy 2005;30:2796-817. 

[67] Klein SA. Calculation of flat-plate loss coefficients. Solar Energy 1975;17: 
79-80. 

[68] Hottel HC, Woertz BB. Performance of flat plate solar heat collectors. Trans- 
actions of the ASME 1942;64:91-104. 

[69] Paisarn N. Effect of porous media on the performance of double-pass flat plate 
solar air heater. International Communication of Heat and Mass Transfer 
2005;32. N 3. 

[70] Sopain K, Supranto, Daud WRW, Othman MY, Yatim B. Renewable Energy 
1999;18:557. 

[71] Paisarn N. On the performance and entropy generation of the double-pass solar 
air heater with longitudinal fins. Renewable Energy 2005;30(9):1345-57. 

[72] Chandra R, Singh NP, Sodha MS. Thermal performance of triple-pass solar air 
collector. Energy Conversion and Management 1990;30(1):41-8. 

[73] Sisson LE, Pitts DR. Elements of transport phenomena. New York: McGraw- 
Hill; 1972. pp. 681-686. 

[74] Ezeike GOI. Energy in Agriculture 1986;5:1-20. 

[75] Choudhury C, Chauhan PM, Garg HP, Garg SN. Cost-benefit ratio of triple pass 
solar air heaters. Energy Conversion and Management 1996;37(1):95-116. 

[76] Ho CD, Yeh HM, Hsieh SM. Improvement in device performance of multi-pass 
flat-plate solar air heaters with external recycle. Renewable Energy 
2005;30:1601-21. 

[77] Pramuang S, Exell RHB. Transient test of a solar air heater with a compound 
parabolic concentrator. Renewable Energy 2005;30:715-28. 

[78] Tchinda R. Thermal behaviour of solar air heater with compound parabolic 
concentrator. Energy Conversion and Management 2007. doi: 10.1016/ 
j.enconman.2007.08.004. 

[79] Hsieh CK. Thermal analysis of CPC collectors. Solar Energy 1981;27:19-27. 

[80] Rabl A, Goodman NB, Winston R. Practical design considerations for CPC solar 
collectors. Solar Energy 1979;22:373-81. 

[81] Brodkey RS, Hershey HC. Transport phenomena. McGraw-Hill Book Com- 
pany; 1988. 

[82] Karim Azharul Md, Hawlader MNA. Development of solar air collectors for 
drying applications. Energy Conversion and Management 2004;45:329-44. 

[83] Azharul Md K, Hawlader MNA. Performance evaluation of v-groove solar air 
collector for drying applications. Applied Thermal Engineering 2006;26: 
121-30. 

[84] Bashria AA, Adam NM, SAPUAN SM, Daud M, Omar H, Megat HM, Abas F. 
Prediction of thermal performance of solar air heaters by internet-based 
simulation. Proceedings of the Institution of Mechanical Engineers Part A 
Journal of Power and Energy 2004;218:579-87. 

[85] Bashria AA, Yousef, Adam NM, Daud M, Omar H, Ahmed MHM, Musa Elradi A. 
Expert system based prediction of the performance of solar air collectors in 
Malaysia. JEE 2004;3:91-101. 

[86] Bashria AA, Yousef, Adam NM. Performance analysis for v-groove absorber. 
Suramaree Journal of Science and Technology 2007;14(1):39-52. 

[87] Verma R, Chandra R, Garg HP. Technical note-optimization of solar air heaters 
of different designs. Renewable Energy 1992;2(4/5):521-31. 

[88] Choudhury C, Garg HP. Performance of air-heating collectors with packed 
airflow passage. Solar Energy 1993;30(3):205-21. 

[89] Paul B, Saini JS. Optimization of bed parameters for packed bed solar energy 
collection system. Renewable Energy 2004;29:1863-76. 

[90] Dixon AG, Cresswell DL. Theoretical prediction of effective heat transfer 
parameters in packed beds. American Institute of Chemical Engineers 
1979;25:663. 

[91] Varshney L, Saini JS. Heat transfer and friction factor correlations for rec- 
tangular solar air heater duct with wire mesh screen matrices. Solar Energy 
1998;62(4):255-62. 

[92] Collier RK. The characterization of crushed glass as a transpired air heating 
solar collector material. In: Proceedings of the International Solar Energy 
Society; 1979. p. 264-8. 

[93] Choudhury C, Garg HP. Performance calculations for closed-loop ai-t-heater 
solar hybrid heating systems with and without a rock bed in the solar air 
heater. Renewable Energy 1993;3(8):897-905. 


R. Tchinda/ Renewable and Sustainable Energy Reviews 13 (2009) 1734-1759 1759 


[94] Choudhury C, Garg HP. Thermal performance of a solar hybrid domestic hot 
water system. Energy 1992;17(7):703-11. 

[95] Duffie JA, Beckmann WA. Solar energy thermal processes. New York: Wiley; 
1974. 


[96] Ramadan MRI, El-Sebaii AA, Aboul-Enein S, El-Bialy E. Thermal performance 
of a packed bed double-pass solar air heater. Energy 2007;32:1524-35. 

[97] Garg HP. Advances in solar energy technology. Dordrecht: Reidel publishing 
Company; 1987. 

[98] Sharma VK, Sharma Sanjay, Mahajan RB, Garg HP. Evaluation of a matrix solar 
air heater. Energy Conversion and Management 1990;30(1):1-8. 

[99] El-Sebaii AA, Abou-Enein S, Ramadan MRI, El-Bialy E. Year round perfor- 
mance of double pass solar air heater with packed bed. Energy Conversion 
and management 2007. 


[100] Jain D, Jain RK. Performance evaluation of an inclined multi-pass solar heater 
with in-built thermal storage on deep-bed drying application. Journal of Food 
Engineering 2004;65:497-509. 

[101] Jain D. Modelling the system performance of multi-tray crop drying using an 
inclined multi-pass solar air heater with in-built thermal storage. Journal of 
Food Engineering 2005;71:44-54. 

[102] Whillier A. Design factors influencing solar collectors, low temperature 
engineering application of solar energy. New York: American Society of 
Heating, Refrigerating and Air conditioning Engineers, Inc.; 1967. 

[103] Lui BYH, Jordan RC. Daily insulation on surfaces tilted towards equator. 
ASHRAE Journal 1962;3(10):53. 

[104] Sharma VK, Rizzi G, Garg HP. Design and development of a matrix type solar 
air heater. Energy Conversion and Management 1991;31(4):379-88. 


